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Abstract 

We review the paradigm of quintom cosmology. This scenario is motivated by the 
observational indications that the equation of state of dark energy across the cos- 
mological constant boundary is mildly favored, although the data are still far from 
being conclusive. As a theoretical setup we introduce a no-go theorem existing in 
quintom cosmology, and based on it we discuss the conditions for the equation of 
state of dark energy realizing the quintom scenario. The simplest quintom model 
can be achieved by introducing two scalar fields with one being quintessence and 
the other phantom. Based on the double-field quintom model we perform a detailed 
analysis of dark energy perturbations and we discuss their effects on current ob- 
servations. This type of scenarios usually suffer from a manifest problem due to 
the existence of a ghost degree of freedom, and thus we review various alternative 
realizations of the quintom paradigm. The developments in particle physics and 
string theory provide potential clues indicating that a quintom scenario may be 
obtained from scalar systems with higher derivative terms, as well as from non- 
scalar systems. Additionally, we construct a quintom realization in the framework 
of braneworld cosmology, where the cosmic acceleration and the phantom divide 
crossing result from the combined effects of the field evolution on the brane and 
the competition between four and five dimensional gravity. Finally, we study the 
outsets and fates of a universe in quintom cosmology. In a scenario with null energy 
condition violation one may obtain a bouncing solution at early times and therefore 
avoid the Big Bang singularity. Furthermore, if this occurs periodically, we obtain 
a realization of an oscillating universe. Lastly, we comment on several open issues 
in quintom cosmology and their connection to future investigations. 
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1 Introduction 



Accompanied by the recent developments on distance detection techniques, such as balloons, 
telescopes and satellites, our knowledge about cosmology has been greatly enriched. These new- 
discoveries in astrophysical experiments have brought many challenges for the current theory of 
cosmology. The most distinguished event is that two independent observational signals on distant 
Type la supernovae (SNIa) in 1998 have revealed the speeding up of our universe [Tf2] . This 
acceleration implies that if the theory of Einstein's gravity is reliable on cosmological scales, then 
our universe is dominated by a mysterious form of matter. This unknown component possesses 
some remarkable features, for instance it is not clustered on large length scales and its pressure 
must be negative in order to be able to drive the current acceleration of the universe. This 
matter content is called "dark energy" (DE). Observations show that the energy density of 
DE occupies about 70% of today's universe. However, at early cosmological epochs DE could 
not have dominated since it would have destroyed the formation of the observed large scale 
structure. These features have significantly challenged our thoughts about Nature. People begin 
to ask questions like: What is the constitution of DE? Why it dominates the evolution of our 
universe today? What is the relation among DE, dark matter and particle physics, which is 
successfully constructed? 



The simplest solution to the above questions is a cosmological constant A [3.4.5,6]. As required 
by observations, the energy density of this constant has to be p\ ~ (10 _3 eV A ) 4 , which seems 
un-physically small comparing to other physical constants in Einstein's gravity. At the classical 
level this value does not suffer from any problems and we can measure it with progressively 
higher accuracy by accumulated observational data. However, questioning about the origin of a 
cosmological constant, given by the energy stored in the vacuum, does not lead to a reasonable 
answer. Since in particle physics the vacuum-energy is associated with phase transitions and 
symmetry breaking, the vacuum of quantum electrodynamics for instance implies a pa about 120 
orders of magnitude larger than what has been observed. This is the worst fine-tuning problem 
of physics. 



Since the fundamental theory of nature that could explain the microscopic physics of DE is 
unknown at present, phenomenologists take delight in constructing various models based on its 
macroscopic behavior. There have been a number of review articles on theoretical developments 
and phenomenological studies of dark energy and acceleration and here we would like to refer 
to Refs. [Tf8"f9"fl0|lll|ll2|ll3|ll4] as the background for the current paper. Note that, the most 
powerful quantity of DE is its equation of state (EoS) effectively defined as wde = Pde/pde, 
where pde and pde are the pressure and energy density respectively. If we restrict ourselves in 
four dimensional Einstein's gravity, nearly all DE models can be classified by the behaviors of 
equations of state as following: 



• Cosmological constant: its EoS is exactly equal to Wx — — 1- 

• Quintessence: its EoS remains above the cosmological constant boundary, that is wq > — 1 
[T5in6] . 

• Phantom: its EoS lies below the cosmological constant boundary, that is wp < —1 [TTfTB"] . 

• Quintom: its EoS is able to evolve across the cosmological constant boundary [TP"] . 
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With the accumulated observational data, such as SNIa, Wilkinson Microwave Anisotropy Probe 
observations (WMAP), Sloan Digital Sky Survey (SDSS) and forthcoming Planck etc., it becomes 
possible in the recent and coming years to probe the dynamics of DE by using parameterizations of 
its EoS, constraining the corresponding models. Although the recent data-fits show a remarkable 
agreement with the cosmological constant and the general belief is that the data are far from being 
conclusive, it is worth noting that some data analyses suggest the cosmological constant boundary 
(or phantom divide) is crossed [T9|20] . which corresponds to a class of dynamical models with 
EoS across —1, dubbed quintom. This potential experimental signature introduced an additional 
big challenge to theoretical cosmology. As far as we know all consistent theories in physics satisfy 
the so-called Null Energy Condition (NEC), which requires the EoS of normal matter not to be 
smaller than the cosmological constant boundary, otherwise the theory might be unstable and 
unbounded. Therefore, the construction of the quintom paradigm is a very hard task theoretically. 
As first pointed out in Ref. [19], and later proven in Ref. [21] (see also Refs. [22, 23, 24, 25 126] ). for 
a single fluid or a single scalar field with a generic lagrangian of form £(4>, (9^0(9^0) there exists 
in general a no-go theorem forbidding the EoS crossing over the cosmological constant. Hence, 
at this level, a quintom scenario of DE is designed to enlighten the nature of NEC violation. Due 
to this unique feature, quintom cosmology differs from any other paradigm in the determination 
of the cosmological evolution and the fate of the universe. 

This review is primarily intended to present the current status of research on quintom cosmology, 
including theoretical constructions of quintom models, its perturbation theory and predictions on 
observations. Moreover, we include the discussions about quintom cosmology and NEC, in order 
to make the nature of DE more transparent. Finally, we examine the application of quintom in 
the early universe, which leads to a nonsingular bouncing solution. 

This work is organized as follows. In Section [2] we begin with the basics of Friedmann- Robertson- 
Walker cosmology and we introduce a concordant model of ACDM, referring briefly to scenarios 
beyond ACDM. In Section [3] we present the theoretical setup of quintom cosmology and we 
discuss the conditions for the DE EoS crossing —1. In SectionH]we introduce the simplest quintom 
scenario, which involves two scalar fields, and we extract its basic properties. Section is devoted 
to the discussion of the perturbation theory in quintom cosmology and to the examination of its 
potential signatures on cosmological observations. Due to the existence in quintom cosmology 
of a degree-of-freedom violating NEC, the aforementioned simplest model often suffers from a 
quantum instability inherited from phantom behavior. Therefore, in Section [6] we extend to a 
class of quintom models involving higher derivative terms, since these constructions might be 
inspired by fundamental theories such as string theory. In Section [7] we present the constructions 
of quintom behavior in non-scalar models, such are cosmological systems driven by a spinor or 
vector field, while in Section [8] we turn to the discussion of quintom-realizations in modified (or 
extended) Einstein's gravity. In Section [9] we examine the violation of NEC in quintom cosmology. 
Additionally, we apply it to the early universe, obtaining a nonsingular bouncing solution in four 
dimensional Einstein's gravity, and we further give an example of an exactly cyclic solution in 
quintom cosmology which is completely free of spacetime singularities. Finally, in Section [10] we 
conclude, summarize and outline future prospects of quintom cosmology. Lastly, we finish our 
work by addressing some unsettled issues of quintom cosmology. Throughout the review we use 
the normalization of natural units c = h = 1 and we define k 2 = 8ttG = M~ 2 . 
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2 Basic cosmology 



Modern cosmology is based on the assumptions of large-scale homogeneity and isotropy of the 
universe, associated with the assumption of Einstein's general relativity validity on cosmological 
scales. In this section we briefly report on current observational status of our universe and in 
particular of DE, and we review the elements of FRW cosmology. 



2.1 Observational evidence for DE 



With the accumulation of observational data from Cosmic Microwave Background measurements 
(CMB), Large Scale Structure surveys (LSS) and Supernovae observations, and the improvement 
of the data quality, the cosmological observations play a crucial role in our understanding of the 
universe. In this subsection, we briefly review the observational evidence for DE. 



2.1.1 Supernovae la 



In 1998 two groups [Tf2] independently discovered the accelerating expansion of our current 
universe, based on the analysis of SNIa observations of the redshift-distance relations. The lumi- 
nosity distance d^, which is very important in astronomy, is related to the apparent magnitude 
m of the source with an absolute magnitude M through 



/i = m — M = 5 log 



di 



10 



Mpc 



+ 25 



(2.1) 



where \i is the distance module. In flat FRW universe the luminosity distance is given by: 



di 



1 + z r z dz' 
H Q Jo E(z') ' 



(2.2) 



with 



E(z) 



H{z) 
Ho 



O r0 (l + zf + fi A0 exp 



l + w(z') 
1 + z' 



1 1 V2 



dz' 



(2.3) 



where Hq is the Hubble constant, w is the EoS of DE component, z is the redshift, and Qio = 
SiTGpio/(3HQ) is the density parameter of each component (with Q m0 + Q r0 + Q\o = 1), i.e its 
energy density divided by the critical density at the present epoch. 



In 1998 both the the high-z supernova team (HSST) [TJ and the supernova cosmology project 
(SCP) [2], had found that distant SNIa are dimmer than they would be in dark-matter dominated, 
decelerating universe. When assuming the flat FRW cosmology (fl m o + ^ao — 1)> Perlmutter et 
al. found that at present the energy density of dark-matter component is: 

ftjjg = 0.28 t° °ol O statistical) (identified systematics) , (2.4) 
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Fig. 1. (Color online) Hubble diagram of SNIa measured by SCP and HSST groups. In the bottom graph, 
magnitudes relative to a universe with (T2 m o = 0.3,f2Ao = are shown. From Ref. \27f . 

based on the analysis of 42 SNIa at redshifts between 0.18 and 0.83 [2]. In FigJT] we present 
the Hubble diagram of SNIa measured by SCP and HSST groups, together with the theoretical 
curves of three cosmological models. One can observe that the data are inconsistent with a A = 
flat universe and strongly support a nonzero and positive cosmological constant at greater than 
99% confidence. See Ref. [28 29.30 31 32 33J for recent progress on the SNIa data analysis. 



2.1.2 CMBandLSS 

Besides the SNIa observations, the CMB and LSS measurements also provide evidences for a 
dark-energy dominated universe. The recent WMAP data [31] are in good agreement with a 
Gaussian, adiabatic, and scale-invariant primordial spectrum, which are consistent with single- 
field slow-roll inflation predictions. Moreover, the positions and amplitudes of the acoustic peaks 
indicate that the universe is spatially flat, with —0.0179 < Qko < 0.0081 (95%CL). 

Under the assumption of DE EoS being w = — 1, the so called ACDM model, the 5-year WMAP 
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Fig. 2. (Color online) Two-dimensional contours of the vacuum energy density, Q\, and the spatial 
curvature parameter SIk- From Ref. 

data, when they combine with the distance measurements from the SNIa and the Baryon Acoustic 
Oscillations (BAO) in the distribution of galaxies, imply that the DE density parameter is I^ao = 
0.726 ± 0.015 at present. In Figf5]we illustrate the two dimensional contour plots of the vacuum 
energy density fl\ and the spatial curvature parameter Q^- One can see that the flat universe 
without the cosmological constant is obviously ruled out. Furthermore, the LSS data of the 
SDSS [35 36. 37] agree with the WMAP data in favoring the flat universe dominated by the DE 
component. 



2.1.3 The age of the universe 

An additional evidence arises from the comparison of the age of the universe with other age- 
independent estimations of the oldest stellar objects. In flat FRW cosmology the age of the 
universe is given by: 

f tu a i — f°° dz f°° dz 

K = Jo Jo H(l + z) = Jo H (l + z)E(z) ' 1 } 

In the dark-matter dominated flat universe (fix = 0, Q m o = 1), relation (I2.5P leads to a universe- 
age: 

r°° dz 2 . 

tu = / , = 777T • ( 2 - 6 ) 

We have neglected the radiation contribution, setting Q y0 = 0, since the radiation dominated 
period is much shorter than the total age of the universe. Combining with the present constraints 
on the Hubble constant from the Hubble Space Telescope Key project (HST) |38j, h = 0.72±0.08, 
the age of universe is 7.4Gyr < t u < lO.lGyr. 

On the other hand, many groups have independently measured the oldest stellar objects and 
have extracted the corresponding constraints the age of the universe: HGyr < t u < 15Gyr 
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Fig. 3. (Color online) The age of the universe t u (in units of H Q 1 ) versus f2 m o- From Ref 

[39|40f41|H2] . Moreover, the present 5- year WMAP data provide the age-limit of t u = 13.72 ± 
0.12Gyr, when assuming the ACDM model. Thus, one can conclude that the age of a flat universe 
dominated by dark matter is inconsistent with these age limits. 

The above contradiction can be easily resolved by the flat universe with a cosmo logical constant. 
Including the DE component, the age of the universe reads: 



dz 



(2.7) 



AO 



Therefore, the age of universe will increase when becomes large. In Figj3] we depict the 
universe age t u (in units of -f^o" 1 ) versus fi m0 . In this case, t u = 1/H m 13.58 Gyr when f2 m0 = 0.26 
and f^Ao — 0.74. As we observe, the age limits from the oldest stars are safely satisfied. 



2.2 A concordant model: ACDM 



From observations such as the large scale distribution of galaxies and the CMB near-uniform 
temperature, we have deduced that our universe is nearly homogeneous and isotropic. Within 
this assumption we can express the background metric in the FRW form, 



ds< 



dt 2 - a 2 (t) 



dr 2 
1 — kr 2 



+ r 2 d£l\ 



(2.8) 
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where t is the cosmic time, r is the spatial radius coordinate, is the 2-dimensional unit 
sphere volume, and the quantity k characterizes the curvature of 3-dimensional space of which 
k = —1,0, 1 corresponds to open, flat and closed universe respectively. Finally, as usual, a(t) is 
the scale factor, describing the expansion of the universe. 



A concordant model only involves radiation p r , baryon matter pi,, cold dark matter pd m arid a 
cosmological constant A. Thus, this class consists the so-called ACDM models. In the frame of 
standard Einstein's gravity, the background evolution is determined by the Friedmann equation 

tt2 8ttG / \ . . 

H = —z-[Pr + Pb + Pdm + Pk + Pa) , (2.9) 



where if = - is the Hubble parameter, and we have defined the effective energy densities of spatial 
curvature and cosmological constant pk = — % ^ Ga 2 and p^ = respectively. Furthermore, the 
continuity equations for the various matter components write as 

Pi + 3#(l + w f > f = 0, (2.10) 

where the equation-of-state parameters Wi = — are defined as the ratio of pressure to energy 

Pi 

density. In particular, they read w r = ~ for radiation, Wb = for baryon matter, u> rfm = for cold 
dark matter, = — | for spatial curvature, and wa = —1 for cosmological constant. One can 
generalize the EoS of the i-th component as a function of the redshift Wi(z), with the redshift 
given as 1 + z — -. Therefore, the evolutions of the various energy densities are given by 



Pi = PiO 



exp|3jf[l + Wi(z)]dln(l + z)} . (2.11) 



From observations we deduce that A is of the order of the present Hubble parameter H , of which 
the energy density is given by 

Pa ~ 10~ 47 GeV 4 . (2.12) 

This provides a critical energy scale 

M A ~ p\ ~ 10" 3 eV . (2.13) 

As far as we know, this energy scale is far below any cut-off or symmetry-breaking scales in 
quantum field theory. Therefore, if the cosmological constant originates from a vacuum energy 
density in quantum field theory, we need to find another constant to cancel this vacuum energy 
density but leave the rest slightly deviated from vanishing. This is a well-known fine-tuning 
problem in ACDM cosmology (for example see Ref. [13] for a comprehensive discussion). 



Attempts on alleviating the fine-tuning problem of ACDM have been intensively addressed in the 
frame of string theory, namely, Bousso-Polchinski (BP) mechanism [44], Kachru-Kallosh-Linde- 
Trivedi (KKLT) scenario [45], and anthropic selection in multiverse [46J. 
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2.3 Beyond ACDM 



DE scenario constructed by a cosmological constant corresponds to a perfect fluid with EoS wa = 
— 1. Phenomeno logically, one can construct a model of DE with a dynamical component, such as 
the quintessence, phantom, K-essence, or quintom. With accumulated astronomical data of higher 
precision, it becomes possible in recent years to probe the current and even the early behavior of 
DE, by using parameterizations for its EoS, and additionally to constrain its dynamical behavior. 
In particular, the new released 5-year WMAP data have given the most precise probe on the 
CMB radiations so far. The recent data-fits of the combination of 5-year WMAP with other 
cosmological observational data, remarkably show the consistency of the cosmological constant 
scenario. However, it is worth noting that dynamical DE models are still allowed, and especially 
the subclass of dynamical models with EoS across —1. 



Quintessence is regarded as a DE scenario with EoS larger than —1. One can use a canonical 
scalar field to construct such models, which action is written as 



S Q = J d 4 x^g 



(2.14) 



By varying the action with respect to the metric, one gets the energy density and pressure of 
quintessence, 

PQ = ^ 2 + W), PQ = \tf-V{4>), (2.15) 
and correspondingly the EoS is given by 



-2V 
+ 2V 



(2.16) 



However, recent observation data indicate that the contour of the EoS for DE includes a regime 
with w < — 1 . The simplest scenario extending into this regime uses a scalar field with a negative 
kinetic term, which is also referred to be a ghost [17J. Its action takes the form 



S, 



and thus the DE EoS is 



Wp 



2 M 



+ 2V 



(2.17) 



(2.18) 



^-2V 

Although this model can realize the EoS to be below the cosmological constant boundary w = —1, 
it suffers from the problem of quantum instability with its energy state unbounded from below 
[50|5"T] . Moreover, if there is no maximal value of its potential, this scenario is even unstable at 
the classical level, which is referred as a Big Rip singularity [18]. 



Another class of dynamical DE model is K-essence 



with its Lagrangian P being a 
general function of the kinetic term [551156] . One can define the kinetic variable X = ^d^^' 
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obtaining 

S K = J d 4 x^P(<P,X) , (2.19) 
and the DE energy density and pressure are given by, 

dP 

PK = 2X—-P, p = Ptf,X). (2.20) 
Consequently, its EoS is expressed as 

w « = - 1 + Wp^p • < 2 ' 21 > 

where the subscript ",X" denotes the derivative with respect to X. By requiring a positive 
energy density, this model can realize either w > — 1 or w < — 1, but cannot provide a consistent 
realization of the — 1-crossing, as will be discussed in the next section. 

Regarding the above three classes of DE scenarios together with the cosmological constant model, 
there is a common question needed to be answered, namely why the universe enters a period of 
cosmic acceleration around today, in which the DE energy density is still comparable to that of 
the rest components. In particular, if its domination was stronger than observed then the cosmic 
acceleration would begin earlier and thus large-scale structures such as galaxies would never have 
the time to form. Only a few dynamical DE models present the so-called tracker behavior [49, 57.9J , 
which avoids this coincidence problem. In these models, DE presents an energy density which 
initially closely tracks the radiation energy density until matter-radiation equality, it then tracks 
the dark matter density until recently, and finally it behaves as the observed DE. 



2.4 Observational Evidence for Quintom DE Scenario 



In this subsection we briefly review the recent observational evidence that mildly favor the 
quintom DE scenario. 

In 2004, with the accumulation of supernovae la data, the time variation of DE EoS was allowed 
to be constrained. In Ref. [20] the authors produced uncorrelated and nearly model independent 
band power estimates (basing on the principal component analysis [58J) of the EoS of DE and 
its density as a function of redshift, by fitting to the SNIa data. They found marginal (2a) 
evidence for w(z) < — 1 at z < 0.2, which is consistent with other results in the literature 



This result implied that the EoS of DE could vary with time. Two type of parameterizations for 
u>de are usually considered. One form (Model A) is: 

wde = u>o + w'z , (2.22) 

where Wq is the EoS at present and w' characterizes the running of wde- However, this parametriza- 
tion is only valid in the low redshift, suffering from a severe divergence problem at high redshift, 
such as the last scattering surface z ~ 1100. Therefore, an alternative form (Model B) was 
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-3.0 -2.5 -2.0 -1.5 -1.0 -0,5 0.0 „3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 



Wo 'WO 

Fig. 4. (Color online) Two-dimensional constraints on the DE parameters in two different parameteri- 
zations. The left panel and right panel correspond to Models A and B respectively (see text). From Ref. 

m 



proposed by Ref. [65|l66] : 

z 

wde = w + wi(l - a) = w + w 1 , (2.23) 

1 ~T Z 

where a is the scale factor and w% = —dw/da. This parametrization exhibits a very good behavior 
at high redshifts. 



In Ref. [19] the authors used the "gold" sample of 157 SNIa, the low limit of cosmic ages and the 
HST prior, as well as the uniform weak prior on Qj^h 2 , to constrain the free parameters of the 
aforementioned two DE parameterizations. As shown in FigJU they found that the data seem to 
favor an evolving DE with EoS being below —1 at present, evolved from w > — 1 in the past. 
The best fit value of the current EoS is wq < —1, with its running being larger than 0. 

Apart from the SNIa data, CMB and LSS data can be also used to study the variation of DE 
EoS. In Ref. [57], the authors used the first year WMAP, SDSS and 2dFGRS data to constrain 
the DE models, and they found that the data evidently favor a strongly time-dependent w^e 
at present, which is consistent with other results in the literature [68 69 ) 7D|71|l72ll73ll74f 75f 76|I77] . 
As observed in Figj5l using the latest 5-year WMAP data, combined with SNIa and BAO data, 
the constraints on the DE parameters of Model B are: w = —1.06 ± 0.14 and w\ = 0.36 ± 0.62 
[341)781)79] . Thus, one deduces that current observational data mildly favor wbe crossing the 
phantom divide during the evolution of universe. However, the ACDM model still fits the data 
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Fig. 5. (Color online) Two-dimensional constraints on the DE parameters (wq,w' ) of Model A (see text). 
From Ref. 

in great agreement. 

3 Quintom cosmology: Theoretical basics 

The scenario that the EoS of DE crosses the cosmological constant boundary is referred as a 
"Quintom" scenario. Its appearance has brought another question, namely why does the universe 
enter a period of cosmic super-acceleration just today. The discussion of this second coincidence 
problem has been carried out extensively in a number of works [19j80 l8T1l82] . However, the explicit 
construction of Quintom scenario is more difficult than other dynamical DE models, due to a 
no-go theorem. 

3.1 A No- Go theorem 

In this subsection, we proceed to the detailed presentation and proof of the "No-Go" theorem 
which forbids the EoS parameter of a single perfect fluid or a single scalar field to cross the 
— 1 boundary. This theorem states that in a DE theory described by a single perfect fluid, or 
a single scalar field minimally coupled to Einstein Gravity with a lagrangian C = P(<f),X), 
in FRW geometry, the DE EoS w cannot cross over the cosmological constant boundary. It has 
been proved or discussed in various approaches in the literature [21,22 23 24,25 26], and in the 
following we would desire to present a proof from the viewpoint of stability of DE perturbations. 
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Let us first consider the fluid-case. Generally, a perfect fluid, without viscosity and conduct heat, 
can be described by parameters such as pressure p, energy density p and entropy S, satisfying a 
general EoS p = p(p, S). According to the fluid properties, such a perfect fluid can be classified 
into two classes, namely barotropic and non-barotropic ones. 



Working in the conformal Newtonian gaugdjj, one can describe the DE perturbations in Fourier 
space as 



S' = -(l + w)(6 -3$') -3H(c 2 s -w)8 (3.1) 

& = -H(l - 3w)6 - -^—6 + k 2 ( + , (3.2) 
1 + w \l + w J v ' 

where the prime denotes the derivative with respect to conformal time defined as rj = J dt/a,H 
is the conformal Hubble parameter, c 2 s = 5p/5p is the sound speed and 



ik j 5T° 

S = Sp/p, 9 = (3.3) 
p + p 

are the relative density and velocity perturbations respectively. Finally, as usual the variables $ 
and \& represent metric perturbations of scalar type. 

If the fluid is barotropic, the iso-pressure surface is identical with the iso-density surface, thus the 
pressure depends only on its density, namely p = p(p). The adiabatic sound speed is determined 
by 



2 _ 21 _ _ W ' ( n A \ 

C a — C s adiabatic . w 01//1 i \ ' W -4 / 

p 3H{l + w) 

From this relation we can see that the sound speed of a single perfect fluid is apparently divergent 
when w crosses —1, which leads to instability in DE perturbation. 

If the fluid is non-barotropic, the pressure depends generally on both its density and entropy, 
p = p(p,S). The simple form of the sound speed given in Eq.f l3.3p is not well-defined. From 
the fundamental definition of the sound speed, that is taking gravitational gauge invariance into 
consideration, we can obtain a more general relationship between the pressure and the energy 
density as follows: 



5p = c 2 s 5p , (3.5) 

with the hat denoting gauge invariance. A gauge- invariant form of density fluctuation can be 
written as 



1 We refer to Ref. [83] for a comprehensive study on cosmological perturbation theory. 
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Sp = S p + 3U(p + p)-^ , (3.6) 



and correspondingly a gauge-invariant perturbation of pressure writes 

5p = 5p + 3nc 2 a (p + p)-^ . (3.7) 
Finally, the gauge-invariant intrinsic entropy perturbation T can be described as 

r = ±( S p - c 2 a 5p) = — (Sp- c 2 a 5p) . (3.8) 
wp wp v ' 

Combining Eqs. fl3.5p -f l3.7p . we obtain the following expression, 



Sp = c s Sp + 



k 2 



, 2 2 mp6(l + w) 3H P 6(l+w) P 6w' 
= c > + c s p p w + -p-. (3.9) 

From Eq.f l3.3p . one can see that 9 will be divergent when w crosses —1, unless it satisfies the 
condition 

9w' = k 2 — . (3.10) 
P 

Substituting the definition of adiabatic sound speed c 2 (relation (13 .4p ) and the condition (I3.10p 
into (13 .7p . we obtain Sp = 0. Since T = ^(<5p — c^5p), it is obvious that due to the divergence 
of c 2 at the crossing point, we have to require 5/3 = to maintain a finite T. Thus, we deduce 
the last possibility, that is 5p = and 5p = 0. From relations f l3.6p .( l3T71) this case requires that 
5p = —c 2 m P e ^ +w ^ a nd 5p = — 3H p 9 ^ +w ^ ; anc i thus 5p = c 2 a 5p. Therefore, it returns to the case of 
adiabatic perturbation, which is divergent as mentioned above. 

In conclusion, from classical stability analysis we demonstrated that there is no possibility for a 
single perfect fluid to realize w crossing —1. For other proofs, see Refs. [23|26] . 

Let us now discuss the case of a general single scalar field. The analysis is an extension of the 
discussion in Ref. [25]. The action of the field is given by 

S = J d 4 xy^P^,X) , (3.11) 

where g is the determinant of the metric g^ v . In order to study the DE EoS we first write down 
its energy- momentum tensor. By definition S g ^S = — J d^x^^-T^dg^y and comparing with the 
fluid definition in the FRW universe one obtains: 
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p = -Ti = C, (3.12) 
p = T ° = 2Xp, x -p, (3.13) 

where "x" stands for Using the formulae above, the EoS is given by 

w = P - = V - = -1+ 2XP > X . (3.14) 

p 2Xp tX - V 2Xp :X - V 

This means that, at the crossing point if:*, Xp t x\t* = 0. Since w needs to cross —1, it is required 
that Xp^x changes sign before and after the crossing point. That is, in the neighborhood of t*, 
(t* - e, t* + e), we have 

Xp tX \ 1 *- e -Xp > x\t>+e<0 ■ (3.15) 

Since X = \(f) 2 is non-negative, relation ( 13. 15ft can be simplified as p t x\t*-e • P,x\t*+t < 0. Thus, 
due to the continuity of perturbation during the crossing epoch, we acquire p t x\t* — 0. 

Let us now consider the perturbations of the field. We calculate the perturbation equation with 
respect to conformal time r\ as 



u" - c 2 s V 2 u - 
where we have defined 



- + 3c 2 (ft' - W 

z 



u = , (3.16) 



u = az^- , z = ^J(f)' 2 \p,x\ ■ (3.17) 
A Fourier expansion of the perturbation function u leads to the dispersion relation: 



u 2 = c 2 k 2 -^-3c 2 s ('H'-'H 2 ) , (3.18) 

with c 2 defined as y^. Stability requires c 2 > 0. Since at the crossing point p t x = and p t x 
changes sign, one can always find a small region where c 2 < 0, unless p t x also becomes zero with 
a similar behavior as p t x- Therefore, the parameter z vanishes at the crossing and thus 7i' — "H 2 
in (I3.18P is finite. Hence, when assuming that the universe is fulfilled by that scalar field it turns 
out to be zero, since at the crossing point we have p t x — 0, z = 0. Consequently, if z" ^ at the 
crossing point the term — will be divergent. Note that even if z" = this conclusion is still valid, 

since in this case = '^- r . But z is a non- negative parameter with z = being its minimum and 

"' (4) 

thus z must vanish, therefore *-r is either divergent or equal to where z is also equal to 
zero as discussed above. Along this way, if we assume that the first (n — l)-th derivatives of z 
with respect to r\ vanish at the crossing point and that z^ ^ 0, we can always use the L 'Hospital 
theorem until we find that ^ = ji"J 2 ) > which will still be divergent. Therefore, the dispersion 
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relation will be divergent at the crossing point as well, and hence the perturbation will also be 
unstable. 



In summary, we have analyzed the most general case of a single scalar field described by a 
lagrangian of the form C = £(0, d^d^cj)) and we have studied different possibilities of w crossing 
the cosmological constant boundary. As we have shown, these cases can either lead to a negative 
effective sound speed (? s) or lead to a divergent dispersion relation which makes the system 
unstable. 



3.2 Conditions of the —1-crossing 



Let us close this section by examining the realization conditions of quintom scenario. Dynamically, 
the necessary condition can be expressed as follows: when the EoS is close to the cosmological 
constant boundary w = — 1 we must have w'\ w= -i ^ 0. Under this condition we additionally 
require both the background and the perturbations to be stable and to cross the boundary 
smoothly. Therefore, we can achieve a quintom scenario by breaking certain constraints appearing 
in the no-go theorem. 

Since we have seen that a single fluid or a single scalar field cannot give rise to quintom, we 
can introduce an additional degree of freedom to realize it. Namely, we can construct a model in 
terms of two scalars with one being quintessence and the other a ghost field. For each component 
separately the EoS does not need to cross the cosmological constant boundary and so their 
classical perturbations are stable. However, the combination of these two components can lead 
to a quintom scenario. An alternative way to introduce an extra degree of freedom is to involve 
higher derivative operators in the action. The realization of quintom scenario has also been 
discussed within models of modified gravity, in which we can define an effective EoS to mimic 
the dynamical behavior of DE observed today. Furthermore, a few attempts have been addressed 
in the frame of string theory, however, models of this type suffers from the problem of relating 
stringy scale with DE. We will describes these quintom models in the following sections. 

Here we should stress again that in a realistic quintom construction one ought to consider the 
perturbation aspects carefully, since conventionally, dangerous instabilities do appear. The con- 
cordance cosmology is based on precise observations, many of which are tightly connected to 
the growth of perturbations, and thus we must ensure their stability. If we start merely with 
parameterizations of the scale factor or EoS to realize a quintom scenario, it will become too 
arbitrary without the considerations of perturbations. On the other hand, if we begin with a 
scenario described by a concrete action which leads to an EoS across w = —1, we can make a 
judgement on the model by considering both its background dynamics and the stability of its 
perturbations. 
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4 The simplest quintom model with double fields 



As we proved in the previous section, a single fluid or scalar field cannot realize a viable quintom 
model in conventional cases. Consequently, one must introduce extra degrees of freedom or in- 
troduce the non-minimal couplings or modify the Einstein gravity. In recent years there has been a 
large amount of research in constructing models with w crossing —1 



4-1 The model 



To begin with, let us construct the simple quintom cosmological paradigm. It requires the simul- 
taneous consideration of two fields, namely one canonical <fi and one phantom a, and the DE is 
attributed to their combination [19]. The action of a universe constituted of a such two fields is 
iHH2]: 



S 



d ii 



^- 2 r - ~«r^«R< 



W) + ^ v d,ad u a - V a (a) + C M 



(4.1) 



where we have set k 2 = 8^G as the gravitational coupling. The term Cm accounts for the (dark) 
matter content of the universe, which energy density pu and pressure pu are connected by the 
EoS Pm — wmPm- Finally, although we could straightforwardly include baryonic matter and 
radiation in the model, for simplicity reasons we neglect them. 

In a flat geometry, the Friedmann equations read [8Tf82] : 



H 2 



H 



K 



PM + Pcf> + Pa 



Pm + Pm + P<p + P<t> + Pa + Pa 



(4.2) 
(4.3) 



The evolution equations for the canonical and the phantom fields are: 



p^ + 3H(p^ + p^ = 0, (4.4) 
p a + 3H( P(7 +p (7 ) = } (4.5) 

where H = a/a is the Hubble parameter. 

In these expressions, p^ and p^ are respectively the pressure and density of the canonical field, 
while p a and p a are the corresponding quantities for the phantom field. They are given by: 
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P^ = ^ 2 + W), P«=^ 2 -W)> (4-6) 

Pa = + , Pa = -\° 2 ~ VM > ( 4 -7) 

where V^,(0), V^(cr) are the potentials for the canonical and phantom field respectively. Therefore, 
we can equivalently write the evolution equations for the two constituents of the quintom model 
in field terms: 

+ 3^0 + ^ = (4.8) 

a + 3Ha-^±=0. (4.9) 
da 

Finally, the equations close by considering the evolution of the matter density: 

p M + 3H(p M +p M )=0. (4.10) 



As we have mentioned, in a double-field quintom model, the DE is attributed to the combination 
of the canonical and phantom fields: 



PDE = P^ + Pa, PDE=P<f>+Pa, (4.11) 

and its EoS is given by 

_Pde Pcf,+Pa u 10 v 

w DE = = ■ • (4.12) 

Pde P<t> + Pa 

Alternatively, we could introduce the "total" energy density p to t = Pm + p<t> + Pa, obtaining: 

p tot + 3H(l + w tot )p tot = 0, (4.13) 

with 



wtot = P(/} + Pa + Pm = + w a Vt a + w m ^m, (4.14) 

p<t> + Pa + Pm 



with the individual EoS parameters defined as 



P* Pa Pm (A , _v 

w 4> = — , w a = — , w M = — , (4.15) 

P<t> Pa Pm 

and 
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are the corresponding individual densities. These constitute: 



^ + Q a = Q DE , (4.17) 

and 

n^ + ^ + viM = 1. (4.18) 

4-2 Phase space analysis 



In order to investigate the properties of the constructed simple quintom model, we proceed to a 
phase-space analysis. To perform such a phase-space and stability analysis of the phantom model 
at hand, we have to transform the aforementioned dynamical system into its autonomous form 
[9|141II142II143] . This will be achieved by introducing the auxiliary variables: 



k4> Kb Ky/VM kJvJo) k J~p^ 

Xa = —f=- , V4, = 7=—— , Va = , z = V- , (4.19) 



together with N = In a. Thus, it is easy to see that for every quantity F we acquire F = Hj^. 
Using these variables we straightforwardly obtain: 



n* = jg = 4 + vl, (4-20) 

^ = j£ = -4 + Vl , (4-21) 
^ ss " 2(P 3 ^ + / ct) = 4 + Vl - + yl , (4.22) 



and 



X 4> y<f> X a Va I A r,o\ 

w <t> = Hri — I ' w ° = ^ — 2 ' ( 4 - 23 ) 

For w iot we acquire: 
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wtot = xj - Vl - xi - y% + (7 - l)z\ (4.25) 

where we have introduced the barotropic form for the matter EoS, defining % = 7 — 1. Finally, 
the Friedmann constraint (14. 2 p becomes: 



xl + yl-x 2 a + yl + z 2 = l. (4.26) 



A final assumption must be made in order to handle the potential derivatives that are present 
in f)4.8p and (14. 9p . The usual assumption in the literature is to assume an exponential potential 
of the form 



V a = V ao e- KX °° , (4.27) 

since exponential potentials are known to be significant in various cosmological models |141|I142|[9] . 
Note that equivalently, but more generally, we could consider potentials satisfying = — 

const and similarly A CT = — ~ cons ^ (f° r example this relation is valid for arbitrary 

but nearly flat potentials |144j ). 

Using the auxiliary variables (T4~T9l . the equations of motion (IP]) . (@~3]), (gSD, ( T4~9l and ( HTTOl) 
can be transformed to an autonomous system containing the variables x^, x a , y,p, y a , z and their 
derivatives with respect to iV = In a. Thus, we obtain: 



X' = f(X), (4.28) 

where X is the column vector constituted by the auxiliary variables, f(X) the corresponding 
column vector of the autonomous equations, and prime denotes derivative with respect to N = 
In a. Then, we can extract its critical points X c satisfying X' = 0. In order to determine the 
stability properties of these critical points, we expand ( I4.28P around X c , setting X = X c + U 
with U the perturbations of the variables considered as a column vector. Thus, for each critical 
point we expand the equations for the perturbations up to the first order as: 



U' = S • U, (4.29) 

where the matrix S contains the coefficients of the perturbation equations. Thus, for each critical 
point, the eigenvalues of H determine its type and stability. 

In particular, the autonomous form of the cosmological system is 
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Label x uc y uc x^c y^ c z c 

A x\ - x 2 a = 1 

DO 1 



Table 1 

The list of the critical points of the simplest quintom model. 



Label rrti m 2 m 3 m 4 Stability 

A -6(l-2) x l c 3(l-^§X a x ac ) 6(l-^)x2 c 3(1 -^V^c) unstable 

B -±(6 + A£) _i(6 + A^) -37 -A^ stable 

C ^ -3(1 - ^) -3(1 - ^) -37 + A| unstable 

D ^ ^ -3(1- |) -3(1 ~i) unstable 



£ ¥ "(3-?) + ^^-.A-^^l unstable 



Table 2 

The eigenvalues and stability of the critical points of the simplest quintom model. 



x; = -3x^ ( 1 + x 2 a - x\ - lz 2 ) + X^y 2 , (4.30) 



\/6 



2 J v 2 



^ = 3^^-4 + ^ + ^ 2 - A 0^^ , ( 4 - 31 ) 
< = -3x CT f 1 + ^ - ^ - Iz 2 ] - A CT ~~~y 2 , (4.32) 



Va = 3ya (^-Xl +Xl + ~Z 2 - \a^X a ) , (4.33) 

z> = 3z(-4 + a* + lz 2 -^) . (4.34) 



2 ; a 2 

A ^ 
2" " Act T 

^ - 1 - 
2 2, 

The real and physically meaningful (i.e corresponding to yi > and < Qi < 1) of them are 
presented in table [U In table [2] we present the eigenvalues of the corresponding matrix H for 
these critical points, which determine their stability properties (a stable point requires negative 
real parts of all eigenvalues). 



From Tables [T]and [21 one can see that the phantom-dominated solution is always a late-time stable 
attractor (stable point B). In particular, this late-time solution corresponds to Qde — &a — 1 an d 
wde = w a = — 1 — A^/3 (see relations (I4.22p and (I4.24p ). i.e to a complete DE domination. This is 
true even if there exists the interaction between the two fields [82J. Thus, if this coupled system is 
initially dominated by the quintessence field, it will eventually evolve into the phantom-dominated 
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phase and the crossing through the phantom divide is inevitable. This behavior provides a natural 
realization of the quintom scenario, and was the motive of the present cosmological paradigm. 
Finally, we mention that during its motion to the late-time attractor, the dynamical system could 
present an oscillating behavior for w E e, with observationally testable effects [52] . 



4-3 State-finder diagnosis 



In this paragraph we are going to present a way to discriminate between the various quintom sce- 
narios, in a model independent manner, following [145] and based in the seminal works [14b'fl47j . 
In those works, the authors proposed a cosmological diagnostic pair {r, s} called statefmder, 
which are defined as 



to differentiate between different forms of DE. In these expressions, q is the deceleration parameter 
q = —ad /a 2 = —d/aH 2 , r forms the next step in the hierarchy of geometrical cosmological 
parameters beyond H and q, and s is a linear combination of r and q. Apparently, the statefmder 
parameters depend only on a and its derivatives, and thus it is a geometrical diagnostic. Since 
different quintom cosmological models exhibit qualitatively different evolution trajectories in the 
s — r plane, this statefmder diagnostic can differentiate between them. 

However, let us use an alternative form of statefmder parameters defined as [145J: 

t . 9 \Ptot + Ptot ;) Ptot [Ptot + Ptot ;) Ptot , A nn \ 

r = 1 + ; — , s = ; — . (4.3b) 

2 Ptot Ptot Ptot Ptot 

Here p to t = Pm + P<t> + Pa + p ra d is the total energy density (note that in order to be more general 
we have also added a radiation part) and p to t is the corresponding total pressure in the universe. 
Since the total energy, the DE (quintom) energy and radiation are conserved separately, we have 
ptot = -3H(p tot +Ptot), Pde = -3H(1 + w D e)pde and p rad = -4Hp rad , respectively. Thus, we 
can obtain 

3 

r = 1 - -[w D e/H - 3w DE (l + w de )}Qde + 2Q md , (4.37) 

-3[w DE /H - 3w DE (l + w DE )]tt DE + 4Q rad 

j (4.38) 



9w DE n DE + 3{l rad 



and 



q = i(l + 3w DE Q DE + Q rad ) , (4.39) 
where as usual fl DE = p DE / ' p tot and Q rad = p rad /ptot- 
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line, dashed line and dot-dashed line respectively), and po the present energy density of our universe. 
The curves r(s) evolve in the time interval ^ 6 [0.5,4] where to *s the present time. Dots locate the 
current values of the statefinder parameters. From Ref. 

We discuss the statefinder for the quintom scenario, imposing various potentials V(<p, a) = 
V^(0) + V a (a). Firstly we assume that there is no direct coupling between the phantom scalar 
field and the normal scalar field, that is we assume the simple exponential potentials of the 
previous subsection: V(a, <fi) = V ao e~ aa + V^e - ^, where a and (3 are constants. As we have 
shown, in this case the universe is moving towards the phantom dominated, late time attractor 
In FigJH] it is depicted the time evolution of statefinder pair {r, s} in the time 
G [0.5,4] where t is the present time [145]. As we can see, in the past and future the 
r — s curve is almost linear, which means that the deceleration parameter changes from one con- 
stant to another nearly with the increasing of time, and the parameters will pass the fixed point 
of LCDM in the future. These trajectories of r(s) are different from other DE models discussed 
in 1147114911501151115211531. 



interval j- 



Similarly, in FigJ7]we can see the r(s) curves for the potentials V(a,4>) = V ao e~ oirj + V^e - ^ + 
y Qe -K.(a+<f>) anc j y^p-^fy = V ao e~ aa + V^e~^ , which lead to late-time attractors. The former 
potential leads to a Big Rip attractor while the latter one to a de Sitter attractor [8182]. Ap- 
parently the upper figure is very similar to Figj6) This shows that for uncoupling and coupling 
exponential potentials the evolutions of our universe are very similar in the time interval we 
consider here. The lower figure is different from Figj6] but has a common characteristic with 
the phantom with power- law potential |148j . quintessence with inverse power-law potential and 
Chaplygin gas model [149] that it reaches the point of LCDM with the increasing of time. This 
is due to the fact that they all lead to the same fate of the universe-de Sitter expansion, but the 
trajectories to LCDM are different, therefore they can be differentiated. 



Let us investigate the case of linear coupling potential V(a, <fi) = k(u + 0) + Acr0, where k and A 
are two constants. The scalar field with a linear potential was first studied in |154j and it has been 
argued that such a potential is favored by anthropic principle considerations [155fl 56 157] and 
can solve the coincidence problem |158j . In addition, if the universe is dominated by quintessence 
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Fig. 7. In the upper figure curves r(s) evolve in the time interval — G [0.5,4] where to is the present 
time for the potential V(o,4>) = V aQ e~ aa + V^e - ^^ + Vbe -K ( CT+< ^. The model parameters are chosen 
as V ao = 0.3/9o, y<j> = 0.6/Jo, Vo = 0.3po, ot = 1, /3 = 1 and k = 1, 2(solid line and dot-dashed line 
respectively). In the lower figure the potential is V(a,(p) = V ao e~ aa2 + V^e - ^^ 1 . The model parameter 
are chosen as V ao = 0.3po, = 0.6po, a = 1 and = 1, 3(solid line and dot-dashed line respectively). 
Dots locate the current values of the statefinder parameters. From Ref. [H5]. 

(phantom) with this potential, it ends with a Big Crunch (Big Rip) [1~8|I159] . The time evolutions 
of statefinder pair {r, s} in the time interval ^ G [0.5,4] are shown in FigJHl We see that in 
the case of negative coupling the diagram is very similar to Figj6] and the upper graph in FigJTl 
which shows that in these cases the evolutions of our universe are similar in the time interval we 
consider here. 

In conclusion, the statefinder diagnostic can differentiate the quintom model with other DE 
models, but it is not very helpful in order to differentiate quintom DE models with some different 
kinds of potentials which lead to a similar evolution of our unverse in the time interval we consider. 

4-4 Cosmic duality 

As we have seen in paragraph 14.11 the simplest quintom model consists of two scalar fields, one 
canonical, quintessence- like and one phantom one [%T|82"] . In |104] it was shown that there exists 
two basic categories of such quintom models: one is quintom-A type, where the canonical field 
dominates at early times while the phantom one dominates at late times and thus the phantom 
divide crossing is realized from above to below, and the other is quintom-B type for which the 
EoS is arranged to change from below —1 to above —1. As we have analyzed in paragraph 14.11 
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Fig. 8. The r — s diagrams for the linear potential case V(a,4>) = k(cj + (p) + Xa<p. The curves r(s) 
evolve in the time interval S [0.5,4]. In the upper graph k = 0.8 and A = —0.2, — 0.8(solid line 
and dashed line respectively). In the lower graph k = 0.8 and A = 0.2, 0.8 (solid line and dashed line 
respectively). Dots locate the current values of the statefinder parameters. From Ref. \145^ - 

the exponential, as well as other simple potential forms, belong to quintom-A type, which is 
realized easily. However, if one desires to construct a specific model of quintom-B type then he 
has to use more sophisticated or fine-tuned potentials, or to add more degrees of freedom |104j . 
Alternatively he can include higher derivative terms (see for example [87]), or add suitably con- 
structed interactive terms which lead to a transition from phantom-to-quintessence domination 
(see [82J). Generally speaking, both quintom-A and quintom-B types could be consistent with 
current observational data. 

One question arises naturally: is there a cosmic duality between these two quintom types? Dual- 
ities in field and string theory have been widely studied, predicting many interesting phenomena 
[160J. The authors of [16111162] have considered a possible transformation with the Hubble pa- 
rameter and they have studied the relevant issues with the cosmic duality |163|I164] . Specifically, 
in [165] a link between standard cosmology with quintessence matter and contracting cosmology 
with phantom, has been shown. Later on this duality was generalized into more complicated DE 
models, where it has been shown to exist too [166 |167lll68|ll69|ll7Ulll71j . In |172|I173] the authors 
have studied this cosmic duality and its connection to the fates of the universe, while in [173] the 
author has also discussed the possibility of realizing the aforementioned duality in braneworld 
cosmological paradigm. A common feature of these studies is that the EoS parameter does not 
cross —1. Therefore, it would be interesting to study the implications of this cosmic duality in 
quintom models of DE, and in particular between quintom-A and quintom-B types. 



We consider the simple quintom model constructed in 14.11 Following 




we can construct 
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a form-invariant transformation by defining a group of new quantities H, p, p and w which keep 
Einstein equations invariant: 



P = p(Pde) , (4.40) 
& = -(-!— Y H . (4.41) 

Under this transformation, we obtain the corresponding changes for the pressure p EE arid the 
EoS w DE , 



p= -?-W)' {pde+pde) iL> {aa2) 



From relations (I4.42p and (14.431) one can see that for a positive , one would be able to 
establish a connection between the quintom-A and quintom-B types. Assuming without loss of 
generality, and as an example for a detailed discussion, that p = p EE in ( 14.421) and (14.431) . we 
can obtain the dual transformation: 



H = -H, (4.44) 

p = -2p DE -p DE , (4.45) 

w = -2- w DE . (4.46) 

Consequently, using the canonical and phantom energy density definitions (14. 6 j) and ( 14. 7p . we 
can extract the dual form of the (quintom) DE Lagrangian: 



£ = \d^a - \d^d^ - Sd{<f>) - 5C 2 {a) , (4.47) 
where 8C\ and 5C 2 are 



5£ 1 = V^) + <P 2 , (4.48) 
5C 2 = V a (a)-& 2 . (4.49) 

Therefore, we can easily see that if the original Lagrangian is a quintom-A type then the dual one 
is a quintom-B type, and vice versa. Thus, under this duality one expects a general connection 
amongst different fates of the universe, and it might be possible that the early universe is linked 
to its subsequent epochs. 

For a specific discussion let us impose a special form for the potentials [105J: 



25 



V^<p) oc -3V2(f) + 2e-^ , (4.50) 
V^a) oc^a 2 + Aa . (4.51) 

Solving explicitly the cosmological equation (14. 2p . (14. 3p . (14. 8p . (14. 9p . neglecting the matter density, 
we study the two periods of the universe evolution. For early times where \t\ <C 1 (in Planck 
mass units), we choose the initial conditions by fixing — > — oo and a — > 0. With these initial 
conditions we can see that the dominant component in DE density is the exponential term of 
V^(0) in (14.501) . and thus the contribution from the phantom potential in (I4.5ip and the linear 
part of quintessence potential in (I4.50p can be neglected. Therefore, the universe behaves like 
being dominated by the quintessence component and it evolves following the approximate 
analytical solution |105j : 



<p~V2\n\t\, a~-t 2 , H~-. (4.52) 

Thus, we see that the scale factor in this period would variate with respect to time as: a oc ±t, 
in which the signal is determined by the positive definite form of the scale factor. Therefore, 
the scale factor here corresponds to the Big Bang or Big Crunch of quintessence-dominanated 
universe. 



The dual form of the solution above is a description of a universe dominated by a phantom 
component with a Lagrangian given by (I4.47P and 



5d + 5C 2 = 2 - o 2 + V^{4>) + V a (a) oc (-3^20 + 4e -V5 * + 2a + ^a 2 ^j . (4.53) 

In addition, the dual Hubble parameter is of the form H ~ — | and for the scale factor we 
acquire a oc ±-. Accordingly, the scale factor of the dual form is tending towards infinity in the 
beginning or the end of universe. From what we have investigated so far, we can see that, for 
the positive branch there is a duality between an expanding universe with initial singularity at 
t = + and a contracting one that begins with an infinite scale factor at t = + . However, for 
the negative branch there is a duality between a contracting universe ending in a big crunch at 
t = 0~ and an expanding one that ends in a final Big Rip at t = 0~. The latter is dominated by a 
phantom component. Besides, in general, under phantom domination the contracting solution is 
not stable, because the phantom universe will evolve into Big Rip or Big Sudden or will expand 
forever approaching to a de Sitter solution, and thus it will not be able to stay in the contracting 
phase forever |175|176j . This problem, however, can be avoided in quintom cosmology since in 
the dual universe with quintom-B DE the increase of kinetic energy of phantom during the 
contraction can be set off by that of quintessence at late time. 

For times where \t\ 3> 1, the phantom component in the first quintom model will dominate and 
the universe will expand. For the specific potentials (14 .501) and ( 14 .511) the mass term in V a (a) 
will gradually play an important role in the evolution of DE. Proceeding as above we obtain: 
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~ y/2\n \t\ , <r~ V2t , H ~ t, 



(4.54) 



where we note that the scale factor is a oc exp(y). Consequently, the scale factor would increase 
towards infinity rapidly for the positive branch, while it would start from infinity for the negative 
branch. In this case, the transformed Lagrangian is (14.471) with 

Sd + 5C 2 oc (-3V2<f) + 4e"^ + 4a + ^a 2 - 2^ . (4.55) 

The universe is evolving with a Hubble parameter H ~ —t and a scale factor a oc exp(— y), which 
is close to singularity related to the origin and the fate of universe. Finally, the component which 
dominates the evolution of the universe is a, that is the scalar field resembling quintessence in 
(14.471) . Consequently we conclude that, for the positive branch, there is a dual relation between an 
expanding universe with a fate of expanding for ever with t — > +oo and a contracting universe 
with a destiny of shrinking for ever with t — > +oo. Meanwhile, for the negative one, there is 
a duality connecting a contracting universe starting from near infinity with t — >■ -co and an 
expanding universe originating from infinity with t — > — oo. 

Having presented the analytical arguments for the duality between quintom A and quintom 
B types, we proceed to numerical investigation. In Figj9] we depict the evolution of the EoS 
parameter of the quintom model and its dual. One can see from this figure that under the 
framework of the duality studied above, the EoS of the quintom model and its dual are symmetric 
around w = —1. Accordingly, in this case quintom A is dual to quintom B rigorously, which 
supports our analytical arguments above. 
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Fig. 9. (Color Online) Evolution of EoS parameter w of the quintom model and its dual as a function 

of the scale factor In a for V = -3V2M 3 4> + 2M 4 e"T + |M 2 <r 2 + 4M 3 a and M is the Planck mass. 
From Ref JTOh^ . 
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In FigJTU]we assume the potentials V^(0) and V a (a) to be exponentials and one can see that the 
EoS parameter for quintom A approaches to a fixed value which corresponds to the attractor 
solution of this type of model [82]. Through the duality we can see that there exists a corre- 
sponding attractor of the quintom B model dual to the former one. finally, in FigJTT] we provide 
another example for the duality. 
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Fig. 10. (Color Online) Evolution of EoS parameter w of the quintom model and its dual as a function 
of the scale factor In a for V = Vo(e~Ai + e~ m) and M is the planck mass. From Ref. 



In summary, we have shown that the quintom model has its dual partner, specifically the quintom 
A model is dual to the quintom B one, while the cosmological equations are form-invariant. These 
two models describe two different behaviors of the universe evolution with one in the expanding 
phase and the other in the contracting one, depending on the imposed initial conditions. The 
cosmic duality, which connects the two totally different scenarios of universe evolution, preserves 
the energy density of the universe unchanged but it transforms the Hubble parameter. 



5 Perturbation theory and current observational constraints 



In this section we investigate the perturbations of quintom DE scenario and the effects of these 
perturbations on current observations. As proved in Section [21 it is forbidden for a single fluid or 
scalar field to realize a quintom scenario in conventional cases, and thus one is led to add extra 
degrees of freedom. Therefore, it is important to check the consistency of this extension at the 
classical level and in particular to analyze the behavior of perturbations when the EoS crosses 
the cosmological constant boundary [2"5] . 
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Fig. 11. (Color Online) Evolution of EoS parameter w of the quintom model and its dual as a function 
of the scale factor In a for V = mi 2 </> 2 + m^a 2 where m\ corresponds to the quintessence component 
mass and ni2 the phantom component mass. From Ref. [1051. 

5.1 Analysis of perturbations in quintom cosmology 



As we have argued above, the quintom scenario needs extra degrees of freedom to the conventional 
models of a single scalar field, and the simplest realization of the quintom is a model with two 
scalar fields or two "effective" scalar fields in the case of higher-derivative operators. In the 
following discussions on quintom perturbations we will restrict ourselves to the double-field model 
of quintom with the lagrangian 

C = C Q + C P , (5.1) 

where 

£g = ^iS^i-^i(0i) (5.2) 



describes the quintessence component, and 

C P 



(5.3) 



for the phantom component. The background equations of motion for the two scalar fields 
1, 2) are 



+ ± a 2 







(5.4) 



where the positive sign is for the quintessence and the minus sign for the phantom. In general 
there will be couplings between the two scalar fields, but for simplicity we neglect them. 

For a complete perturbation study, both the fluctuations of the fields, as well as those of the 
metric, need to be considered. In the conformal Newtonian gauge the perturbed metric is given 
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by 

ds 2 = a 2 (r)[(l + 2^)dr 2 - (1 - 2§)dx i dx i ] . (5.5) 

Using the notations of [83], the perturbation equations satisfied by each component of the two- 
field quintom model read: 



Si 

0; 

where 



ll + Wi)(0i-3$)-3H 



'SPi 
dpi 



-H(l-3wi)9i- 



Wi 



l + w, 



6i + k 2 



Wij ^ , 
( SPi/Spi 

[1 + Wi 



(5.6) 
(5.7) 



$i = (k 2 /4>i)54>i , o-j = , Wi 



(5.8) 



and 



5P = 6 Pi - 2V!5<t>i = 5 Pl + ^ [37/(1 - w?) + w< 



Combining Eqs. ( 15. 6ft . ( 15. 7ft and ( 15. 9ft . we have 



0i = 2m 



A- 2 



l + w. 



■Si + k 2 ^ 



5 t = -(1 + Wi)(0i - 3$) - 3H(l - Wi)Si - m 



Wi + m{\ - w 2 ) 
k 2 



9i . 



(5.9) 



(5.10) 
(5.11) 



Since the simple two-field quintom model is essentially a combination of a quintessence and 
a phantom field, one obtains the perturbation equations of quintom by combining the above 
equations. The corresponding variables for the quintom system are 



w 



quintom 
S quintom 
6 quintom 



Hi pi 
J2i PiSj 
Eipi 

Hi(Pi+Pi)Qi 
EiiPi + P) 



(5.12) 
(5.13) 
(5.14) 



Note that for the quintessence component, — 1 < w\ < 1, while for the phantom component, 
w 2 < —I. 

The two-field quintom model is characterized by the potentials V*. Let us consider Vi(<pi) = 
^m 2 (j) 2 . In general the perturbations of (pi today stem from two origins, the adiabatic and the 
isocurvature modes. If we use the gauge invariant variable Q = — $ — instead of Si, and the 
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relation $ = \l/ in the universe without anisotropic stress, the equations (15. lip and (I5.10p can be 
rewritten as, 



(5.15) 



6 l = 2W l + k 2 (3Q + 4<S>) , 



(5.16) 



where 



Ci 



1 + Wi 



+ 3H(1 -Wi) = d Q [ln(a 6 \ Pi + Pl \)} . 



(5.17) 



( a is the curvature perturbation on the uniform- density hypersurfaces for the a-component in the 
universe |177j . Usually, the isocurvature perturbations of 4>i are characterized by the differences 
between the curvature perturbation of the uniform-0j-density hypersurfaces and that of the 
uniform-radiation-density hypersurfaces, 



where the subscript r represents radiation. We assume there are no matter isocurvature pertur- 
bations and thus Cm = Cr- Eliminating Q from equations (I5.15P and (I5.16p . we obtain a second 
order equation for (9j, namely 



This is an inhomogeneous differential equation and its general solution is the sum of the general 
solution of its homogeneous part and a special integration. In the following, we will show that the 
special integration corresponds to the adiabatic perturbation. Before the era of DE domination, 
the universe is dominated by some background fluids, for instance, the radiation or the matter. 
The perturbation equations of the background fluid are, 





Bi + (Ci - m% + (CM -2U + k 2 )9i = k 2 



(4$ + Ci$) . 



(5.19) 



Cf = -Of/3 , 




(5.20) 



From the Poisson equation 




(5.21) 



on large scales we approximately acquire: 
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$ ~ -c, - y 2 6 f . (5.22) 
Combining these equations with H = 2/[(l + 3u>/)t], we obtain (note numerically Of ~ 0(k 2 )Q) 



Cf 

Of 



5 + 3wf 
~3(l+w f ) 
k 2 {l + 3w f ) 
3(1 + «;/) 



$ = const. 



(5.23) 



Thus, we can see from Eq. (15.191) that there is a special solution, which is given approximately 
on large scales by 0°/ d = Of, and from Eq. (I5.16P we have Q d = £*/■ This indicates that the 
special integration of ( 15.19)) corresponds to the adiabatic perturbation. Hence, concerning the 
isocurvature perturbations of (pi, we can consider only the solution of the homogeneous part of 
(EH}: 



6\ + (C< - 2H)6i + {C{H -2U + k 2 )0i = . 
These solutions are represented by Of and Q so and the relation between them is 



(5.24) 



Si 



0j 



2Wi 



3k 2 



Since the general solution of Q is = Q d + C 
simply Si r = 3C S °- 



I, so 



(5.25) 



Cr + C SO ; t ne isocurvature perturbations are 



In order to solve Eq. (I5.24p we need to know the forms of Cj and H as functions of time r. 
For this purpose, we solve the background equations (15.41) . In radiation dominated period, a = 
At , H = 1/t and we have 



-1/2 



-1/2 



A A 
AiJi/^—mxT 2 ) + A 2 J~i/a(— vn X T 2 ) 

A ~ A 

A 1 I l/i (—m 2 T 2 ) + A 2 I-i/ 4 (— m 2 T 2 ) 



(5.26) 
(5.27) 



respectively, where A, A, and A\ are constants, J„(x) is the z/th order Bessel function and I v (x) 
is the vt\i order modified Bessel function. Usually the masses are small in comparison with the 
expansion rate in the early universe m« <C rl/a, and we can approximate the (modified) Bessel 
functions as J v (x) ~ x v (c\ + c 2 x 2 ) and I v (x) ~ x p (c\ + c 2 x 2 ). We note that J_i/4 and J_i/4 are 
divergent when x — >■ 0. Given these arguments one can see that this requires large initial values 
of (pi and 4>i if A 2 and A 2 are not vanished. 

If we choose small initial values, which is the natural choice if the DE fields are assumed to survive 
after inflation, only A\ and Ai modes exist, so will be proportional to r 3 in the leading order. 
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Thus, the parameters Cj in equation (15.171) will be Cj = 10/r (we have used \pi + Pi\ = (pf/a 2 ). 
Therefore, we obtain the solution of Eq. (15.241) as 



T 



Dn cos(kr) + D l2 sin(fcr) 



(5.2* 



That is, 9\ so presents an oscillatory behavior, with an amplitude damping with the expansion of 
the universe. The isocurvature perturbations Q so decrease rapidly. 

On the other hand, if we choose large initial values for (pi and (pi, A 2 and A 2 modes are present, 
<pi will be proportional to r~ 2 in the leading order and Cj = 0. Thus, the solution of ( 15. 24ft is 



Of 



T 



Dn cos(A;r) + D i2 sin(fcr) 



(5.29) 



Therefore, 9\ so will oscillate with an increasing amplitude, so Q so remains constant on large scales. 

Similarly, during matter dominated era, a = Bt 2 , 7-L = 2/r, and the solutions for the fields (pi 
read: 



>1=T 



-3 



> 2 = r 



-3 



B B 
Bi sin( — miT 3 ) + B 2 cos( — mir 3 ) 
3 3 



B B 
B\ sinh( — m 2 r 3 ) + B 2 cosh( — m 2 T 3 ) 
3 3 



(5.30) 
(5.31) 



We acquire the same conclusions as those reached by the aforementioned analysis for the radiation 
dominated era. If we choose small initial values at the beginning of the matter domination we 
find that the isocurvature perturbations in (pi decrease with time. On the contrary, for large 
initial values the isocurvature perturbations remain constant at large scales. This conclusion is 
expectable. In the case of large initial velocity, the energy density in the scalar field is dominated 
by the kinetic term and it behaves like the fluid with w = 1. The isocurvature perturbation in 
such a fluid remains constant on large scales. In the opposite case, however, the energy density 
in the scalar field will be dominated by the potential energy due to the slow rolling. It behaves 
like a cosmo logical constant, and there are only tiny isocurvature perturbations in it. 

In summary, we have seen that the isocurvature perturbations in quintessence-like or phantom- 
like field with quadratic potential, decrease or remain constant on large scales depending on 
the initial velocities. In this sense the isocurvature perturbations are stable on large scales. 
The amplitude of these perturbations will be proportional to the value of Hubble rate evaluated 
during the period of inflation H in f, if their quantum origins are from inflation. For a large H in f the 
isocurvature DE perturbations may be non-negligible and will contribute to the observed CMB 
anisotropy |178lll79j . In the cases discussed here, however, these isocurvature perturbations are 
negligible. Firstly, large initial velocities are not possible if these fields survive after inflation as 
mentioned above. Secondly, even if the initial velocities are large at the beginning of the radiation 
domination, they will be reduced to a small value due to the small masses and the damping 
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effect of Hubble expansion. In general the contributions of DE isocurvature perturbations are 
not very large |180j , and here for simplicity we have assumed that H in f is small enough and thus 
the isocurvature contributions are negligible. Therefore, in the next subsection we focus on the 
effects of the adiabatic perturbations of the quintom model with two scalars fields. 



5.2 Signatures of perturbations in quintom DE 



Based on perturbation equations (15.131) and f)5.14p . we modify the code of CAMB [181] and we 
study preliminarily in this subsection the quintom observational signatures. For simplicity we 
impose a flat geometry as a background, although this is not necessary. Moreover, we assume the 
fiducial background parameters to be fib = 0.042, Qdm = 0.231, Qde = 0.727, where b stands for 
baryons, DM for dark matter and DE for dark energy, while today's Hubble constant is fixed at 
H = 69.255 km/s Mpc~ 2 . We will calculate the effects of perturbed quintom on CMB and LSS. 

In the two-field quintom model there are two parameters, namely the quintessence and phantom 
masses. When the quintessence mass is larger than the Hubble parameter, the field starts to 
oscillate and consequently one obtains an oscillating quintom. In the numerical analysis we fix 
the phantom mass to be mp ~ 4 x 10~ 61 M P , and we vary the quintessence mass with the 
typical values being wiq = 2 x 1CT 61 M P and 8 x 10~ 61 M p respectively. In Fig. [12] we depict 
the equation-of-state parameters as a function of the scale factor, for the aforementioned two 
parameter-sets, and additionally their corresponding effects on observations. We clearly observe 
the quintom oscillating behavior as the mass of quintessence component increases. After reaching 
the w — — 1 pivot for several times, w crosses —1 consequently with the phantom-component 
domination in dark energy. As a result, the quintom fields modify the metric perturbations 
Sgoo = 2a 2 5ga = 2a 2 §5ij, and consequently they contribute to the late-time Integrated Sachs- 
Wolfe (ISW) effect. The ISW effect is an integrant of $ + over conformal time and wavenumber 
k. The above two specific quintom models yield quite different evolving $ + \1/ as shown in the 
right graph of FigJT2l where the scale is k ~ 10~ 3 Mpc -1 . As we can see, the late time ISW 
effects differ significantly when DE perturbations are taken into account (solid lines) or not 
(dashed lines). 

ISW effects play an important role on large angular scales of CMB and on the matter power 
spectrum of LSS. For a constant EoS of phantom the authors of [182] have shown that the low 
multipoles of CMB will be significantly enhanced when DE perturbations are neglected. On the 
other hand, for a matter-like scalar field where the EoS is around zero, perturbations will also 
play an important role on the large scales of CMB, as shown in [183]. The results on CMB and 
LSS reflect the two combined effects of phantom and oscillating quintessence. Note that in the 
early studies of quintessence effects on CMB, one usually considers a constant w e ff instead 

_ Jdan(a)w(a) 
We " = Jdan(a) ' (5 ' 32) 

however this is not enough for the study of effects on SNIa, nor for CMB, when the EoS of 
DE has a very large variation with redshift, such as the model of oscillating quintom considered 
above. 
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Fig. 12. (Color online) Effects of the two-field oscillating quintom on the observables. The mass of the 
phantom field is fixed at 4 x 10 _61 M P and the mass of the quintessence field are 2 x 10~ 61 M P (thicker 
line) and 8 x 1CT 61 M P (th inner line) respectively. The upper right graph illustrates the evolution of 
the metric perturbations + \P of the two models, with (solid lines) and without (dashed lines) DE 
perturbations. The scale is k ~ 10~ 3 Mpc^ 1 . The lower left graph shows the CMB effects and the lower 
right graph delineates the effects on the matter power spectrum, with (solid lines) and without (dashed 
lines ) DE perturbations. From Ref. jl25f . 

To analyze the oscillating quintom-model under the current observations, we perform a pre- 
liminary fitting to the first year WMAP TT and the TE temperature-polarization cross-power 
spectrum, as well as the recently released 157 "Gold" SNIa data [28]. Following Refs. |184|I185] in 
all the fittings below we fix r = 0.17, n^/t 2 = 0.135 and Q^h 2 = 0.022, setting the spectral index 
as ns = 0.95, and using the amplitude of the primordial spectrum as a continuous parameter. In 
the fittings of oscillating quintom we have fixed the phantom-mass to be mp ~ 1.2 x 10~ 61 M P . 
FigH2] delineates 3cr WMAP and SNIa constraints on the two-field quintom model, and in ad- 
dition it shows the corresponding best fit values. In the left graph of Fig{13] we present the 
separate WMAP and SNIa constraints. The green (shaded) area marks WMAP constraints on 
models where DE perturbations have been included, while the blue area (contour with solid 



35 




(m / m, : ) (m / m„) 



Fig. 13. (Color online) 3a WMAP and SNIa constraints on two-field quintom model, shown together 
with the best fit values. rriQ and mp denote the quintessence and phantom mass respectively. We have 
fixed mp ~ 1.2 x 10~ 61 M p and we have varied the value of tuq. Left graph: separate WMAP and SNIa 
constraints. The green (shaded) area marks the WMAP constraints on models where DE perturbations 
have been included, while the blue area (contour with solid lines) corresponds to the case where DE 
perturbations have not been taken into account. Right graph: combined WMAP and SNIa constraints 
on the two-field quintom model with perturbations (green/ shaded region) and without perturbations (red 
region/ contour with solid lines). From Ref. ]25^ . 

lines) is the corresponding area without DE perturbations. The perturbations of DE have no 
effects on the geometric constraint of SNIa. The right graph shows the combined WMAP and 
SNIa constraints on the two-field quintom model with perturbations (green/shaded region) and 
without perturbations (red region/contour with solid lines). We conclude that the confidence 
regions indeed present a large difference, if the DE perturbations have been taken into account 
or not. 

So far we have investigated the imprints of oscillating quintom on CMB and LSS. Now we consider 
another example, in which w crosses —1 smoothly without oscillation. It is interesting to study 
the effects of this type of quintom model, with its effective EoS defined in (15.321) exactly equal to 
— 1, on CMB and matter power-spectrum. This investigation will help to distinguish the quintom 
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Fig. 14. (Color online) Comparison of the effects of the two-field quintom model with w e ff = — 1 and 
of the simple cosmological constant, in CMB (WMAP), the metric perturbations + VP (the scale is 
k ~ 10 -3 Mpc^ 1 ) and the linear growth factor. The binned error bars in the upper right graph are 
WMAP TT and TE data ]18(^187j . From Ref. 

model from the cosmological constant. We have realized such a quintom model in the lower right 
graph of Fig. [TU which can be easily given in the two-field model with lighter quintessence mass. 
In this example we have set uiq ~ 5.2 x 1CT 62 M P and mp ~ 1.2 x 10~ 61 M P . Additionally, we 
assume that there is no initial kinetic energy. The initial value of the quintessence component is 
set to (fin = 0.045M p , while for the phantom part we impose 02« — 1-32 x 10~ 3 M p . We find that 
the EOS of quintom crosses — 1 at z ~ 0.15, which is consistent with the latest SNIa results. 

The quintom scenario, which is mildly favored by current SNIa data, needs to be confronted with 
other observations in the framework of concordance cosmology. Since SNIa offer the only direct 
detection of DE, this model is the most promising to be distinguished from the cosmological 
constant and other dynamical DE models which do not get across —1, by future SNIa projects 
on the low redshift (for illustrations see e.g. [20] ). This is also the case for quintom scenario in the 
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full parameter space: it can be most directly tested in low redshift Type la supernova surveys. 
In the upper left graph of Fig. [TD we delineate the different ISW effects among the cosmological 
constant (red/light solid) and the quintom model which exhibits w e ff = — 1, with (blue/dark 
solid) and without (blue dashed) perturbations. Similarly to the previous oscillating case, the 
difference is big when switching off quintom perturbations and much smaller when including the 
perturbations. In the upper right graph we find that the quintom model cannot be distinguished 
from a cosmological constant in light of WMAP. The two models give almost exactly the same 
results in CMB TT and TE power-spectra when including the perturbations. We deduce that 
the difference in CMB is hardly distinguishable even by cosmic variance. 

So far we have seen that CMB observations cannot distinguish between a quintom model with 
w e ff = — 1 and a cosmological constant. Thus, in order to acquire distinctive signatures, we have 
to rely on other observations. To achieve that we need to consider the physical observables which 
can be affected by the evolving w sensitively. In comparison with the cosmological constant, 
such a quintom model exhibits a different evolution of the universe's expansion history, and in 
particular it gives rise to a different epoch of matter-radiation equality. The Hubble expansion 
rate H is 

H = ^ = H [Q M a- 3 + Q r a~ 4 + X] 1 / 2 , (5.33) 



a 



where X, the energy density ratio of DE between the early times and today, is quite different 
between the quintom-QDM. and ACDM. In the ACDM scenario X is simply a constant, while in 
general for DE models with varying energy density or EoS, 

X = Q DE a- 3 e- 3 I w{a)dhia . (5.34) 

The two models will give different Hubble expansion rates. This is also the case between the 
quintom model with w e ff = — 1 in the left graph of Fig. [14] and a cosmological constant. 

Finally, we mention that different H leads directly to different behaviors of the growth factor. In 
particular, according to the linear perturbation theory all Fourier modes of the matter density 
perturbations grow at the same rate, that is the matter density perturbations are independent 
of k: 

5 k + H5 k - AnGa 2 p M 5 k = . (5.35) 

The growth factor -Di(a) characterizes the growth of the matter density perturbations: Di(a) = 
5k{a)/5k(a = 1) and it is normalized to unity today. In the matter-dominated epoch we have 
Di(a) = a. Analytically -Di(a) is often approximated by the Meszaros equation 

5Q M H{a) f a da' 
Dl{a) = —2H^Jo (a>H(a>)/H o) r (5 ' 36) 

Therefore, we can easily observe the difference between the quintom and cosmological constant 



scenarios, due to the different Hubble expansion rates. More strictly one needs to solve Eq. fl5.35p 
numerically. In the lower left graph of Fig. [TH we show the difference of D\ (a) between the 



2 One should notice that the Meszaros equation is an exact solution to the differential equation for the 
linear-theory growth factor only when dark energy is a cosmological constant or it is absent. However, 
it is still a good approximation in the case of dynamical dark energy models on small length scales of 
the universe, as in the analysis of Section 3.1 of Ref. [189J. 
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quintom with w e ff = — 1 and the cosmological constant. The difference in the linear growth 
function is considerably large in the late-time evolution, and thus possibly distinguishable in 
future LSS surveys and in weak gravitational lensing (WGL) observations. WGL has emerged 
with a direct mapping of cosmic structures and it has been recently shown that the method 
of cosmic magnification tomography can be extremely efficient |190|I191|I192] . which leaves a 
promising future for breaking the degeneracy between quintom and a cosmological constant. 



6 Quintom models with higher derivative terms 

As we demonstrated previously, we usually need to introduce new degrees of freedom into a 
normal lagrangian in order to obtain a viable quintom scenario. One approach is to construct a 
double-field quintom model introduced in Section IV. In this section we provide an alternative 
possibility of introducing extra degrees of freedom for the realization of the transition between 
quintessence phase and phantom phase |193j . This model is originally proposed by Lee and Wick 
to address on hierarchy problem in particle physics in the 70's |194|ll95j . and recently applied in 
[87] involving higher derivative operators in the Lagrangian. 



6.1 Lee-Wick model 



In this subsection we present a single field model with the inclusion of higher derivatives, following 
[87] . We consider a Lagrangian of the form 



where X = l/2V At 0V /i and □ = V M V M is the d'Alembertian operator. We mention that this 
model is just an effective theory and one assumes that the operators associated with the higher 
derivatives can be derived from some fundamental theories, for instance due to the quantum 
corrections or the non-local physics in string theory |196|ll97|ll98j . Additionally, with the inclusion 
of higher derivative terms to the Einstein gravity, the theory is shown to be renormalizable 
[199J. Finally, higher derivative operators have been shown to be capable of stabilizing the linear 
fluctuations in the scenario of "ghost condensation" [20U|2Ulf2 02 203 204 |205] . 

We consider the simple model with 



where M is a constant with mass dimension and c is a dimensionless constant. The energy- 
momentum tensor reads [87] : 



C = C 0,X,(D0)(D0),(V / ,V»(V^V^) 



(6.1) 



(6.2) 
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rpflv _ fW 



oVp0V^+^(D0) 2 + — V0V p (D 



M 2 

c 



V 



M 2 



and the equation of motion is given by 



□ c 



M 2 



dV 



. 



However, the energy-momentum tensor (16.31) can be rewritten as 



where x and r/> are defined by 



X = ' ^ = + X • 



M 2 
c 



(6.3) 



(6.4) 



(6.5) 



(6.6) 



It is not difficult to see that the energy-momentum tensor (16.51) can be derived from the following 
Lagrangian 



C 



1 1 M 2 

--Vp^VV + 2 V ^v^x - V{4, - x) - -^x 2 



(6.7) 



with ip and \ being two independent fields. The variations of the Lagrangian (16. 7p . with respect 
to the fields ip and x respectively, give rise to the equations of motions 



_ V' = , 
M 2 

°x + — x - v = o , 

c 

where the prime denotes the derivative with respect to ip — x- 



(6.8) 
(6.9) 



When one imposes ip = <fi + x 111 fl6.7[) and (16.81) . he recovers the equation of motion (16.41) of 
the single field model and the transformation equation ( 16.61) . Thus, one can see that the single 
field model (16. 2p is equivalent to the one with two fields (16. 7ft . where one is the canonical scalar 
field x an d the other is the phantom field if). Therefore, assuming that x dominates the DE 
sector in the early time, the EoS parameter will be larger than —1, but when the phantom 
mode becomes dominant the EoS parameter will become less than —1, crossing the cosmological 
constant boundary at an intermediate redshift. 



In general the potential term V(if) — x) (equivalently V(<p) in our single field model) should 
include interactions between the two fields if) and x- However, for some specific choices of the 
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potential V these two modes could decouple. As an example we consider V = (l/2)m 2 2 . Then 
the Lagrangian (16.71) can be "diagonalized" as: 



1 2i2 "^2/2 



-1 - - V ^2 V^<P2 ~ T^l^l ~ ~ m 



2^2 > 



(6.10) 



through the transformation 



U 



y-a 2 Oi 



(6.11) 



\02 



where 



1 / 4cm 2 \ 1/4 / / 4cm 2 \ 
1 / 4cm 2 \~ 1/4 / / 4cm 2 \ 

a2 = 2( 1+ Mr) (V 1+ m^ + 1 J 

and 



9 M 2 / / 4cm 2 \ 

m ' = ^(V 1 + i^ + 1 J • 

9 M 2 / / 4cm 2 \ 

m ^ = it (v 1 + ip - v • (6 ' 13) 

Hence, one can see that the model (16. 2ft with V = (l/2)m 2 2 is equivalent to the uncoupled 
system (I6.10p . The two modes (pi and 02 evolve independently in the universe. The positivity of 
the parameter c guarantees the positivity of m 2 (m 2 is always positive as long as 1 + ^p- > 0). 
In the limit m <C M, the masses of the two modes are approximately mi ~ Mj y/c and m 2 — m. 
In fact, 0i and 02 are the eigenfunctions of the d'Alembertian operator □, with the eigenvalues 
— m 2 and mj respectively. The solution to the equation of motion (16. 4 j) 0, is decomposed by these 
eigenfunctions as = (1 + ^ff|-)~ 1 ^ 4 (0i +02)- Finally, we mention that the perturbations in the 
phantom component 02 do not present any unphysical instabilities at the classical level. Thus, 
the model is free from the difficulties of the singularity and the gravitational instabilities of the 
general k-essence-like scenarios [87] . 

In Fig. H5] we present the EoS parameter of the single field model as a function of In a, where 
w clearly crosses —1 during evolution. On the same figure it is also depicted the corresponding 
evolution of EoS parameter of the equivalent two-field model. Indeed, one can see that the two 
models coincide. Therefore, the high derivative single scalar field model of DE presents the — 1- 
crossing, and it is equivalent to the usual two-field quintom paradigm. However, note that when 
we consider the interactions between these two fields in the potential, they may possibly show 
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Fig. 15. (Color online) Evolution of the EoS parameters for the high derivative single-field model and 
the equivalent two-field model with one quintessence and one phantom field. From Ref. \87$ . 

some different behaviors. 



One can easily generalize the aforementioned model considering the Lagrangian 



(6.14) 



where x 



^r-D0 and = 4>(ip, x)- Using A((j>) = —1, V((f)) = this Lagrangian is simplified as 



C 

where < 



(6.15) 



— (ip + x)- The equations of motion for the two fields are: 



M 2 v a 2 
Mia n 



Ml 



— — y H -\ 

2C 2 C 







(6.16) 



where the prime is the derivative with respect to ip + X- I n FRW universe, □ = d 2 /dt 2 + 3Hd/dt 
with H being the Hubble expansion rate. Furthermore, one can easily extract the density and 
the pressure of DE as 
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p = -T + T + ^> (6 ' 17) 

p= -T + y-^' (6 - 18) 

and thus the EoS parameter writes 
p 2(x 2 -^ 2 ) 

= ; ^-2-, 6.19 

where the dot represents the derivative with respect to the cosmic time. Again, \ behaves as a 
canonical field while tp as a phantom one. However, in this case they couple to each other through 
an effective potential, that is the last term in the right hand side of equation ( 16. 15(1 : 



Veff = ^H7Z 3i • ( 6 -20) 



M 2 X 2 
2C[^, X )] 

As an example one can consider 



7T 

C[<fi(ip, X)] = Cq<- + arctan 



M v 



(6.21) 



Thus, C is almost constant when \ip + %| 3> 0, and ip and x are nearly decoupled at these 
regimes, which are dubbed "weak coupling" regimes. By contrast, the two fields couple tightly 
in the "strong coupling" regime where \ip + x\ ~ 0. In the weak coupling regime, as shown in 
(I6.16I) . the phantom-like field ip behaves as a massless scalar field and its energy density — (l/2)ip 2 
dilutes as a~ 6 , with a is the scale factor of the universe. The quintessence-like field % has a mass 



term with m x = M v j y C(ip + %). Its behavior is determined by the ratio of m x /H. If m x <C H, x 
is slow- rolling and it behaves like a cosmo logical constant. On the other hand, in cases m x ^> H 
the kinetic term and potential oscillate coherently and their assembly evolves as a -3 , just like 
that of collisionless dust. 

The EoS parameter of such a DE model is depicted in Fig. [161 The initial conditions are chosen 
in order for the DE to start in the strong coupling regime with x > at high red-shift. As we 
observe, DE is frozen quickly at initial times, and then it becomes significant around redshift 
z ~ 1, and the phantom kinetic term ip 2 becomes larger than the quintessence one x 2 - The DE 
crosses the phantom divide from above to below. In the future, the system evolves into the weak 
coupling regime, tp 2 dilutes quickly and x gets a large mass. The whole system behaves as cold 
matter and its EoS oscillate around the point of w — 0. So, the universe will exit from the 
phantom phase and end in the state of matter-domination. This result resembles the late-time 
behavior of "B-inflation" model [206j . Therefore, the present cosmological scenario differs from 
the simple quintom DE model, since it leads to a double crossing of the phantom divide and 
it can give a nice exit from phantom phase, in which it stays just for a finite time, avoiding 
the realization of a Big Rip [T8|129f207ir2 08 : 209j. In addition, at late times w oscillates with 
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Fig. 16. (Color online) The EoS parameter w as a function of 
In a for C = 6.0 x 10 121 , A = 25, and the initial values are 
^ = -0.052M p , 4>i = 1-40 x 10~ 63 M 2 , X i = 0.050M p , Xi = -1-44 x 10" 63 M 2 , with Q DEo « 0.73. 
From Ref. J57p. 

high frequency around zero, so the universe will end in an accelerated expansion at some time. 
Finally, note that the quintom model based on higher derivative theory possesses also ghost 
degrees of freedom, and thus it could face the problem of quantum instability [5Tf50] . However, 
as pointed out in |210j . the problem of quantum instability arises because <fi and □</> are quantized 
in canonical way independently. In fact, in the present model both of them are determined by <p, 
and a more appropriate quantization method seems to be possible to avoid the instability [97J. 



6.2 Quintom DE inspired by string theory 



In recent years, there have been active studies on relating quintom scenario with fundamental 
physics, and in particular with string theory. In such a context a rolling tachyon field in the world 
volume theory of the open string, stretched between a D-brane and an anti-D-brane or a non-BPS 
D-brane, plays the role of the scalar field in the quintom model |101f211|l212ll213f214f215f216f217j . 
The effective action used in the study of tachyon cosmology consists of the standard Einstein- 
Hilbert action plus an effective action for the tachyon field on unstable D-brane or (D-brane)- 
(anti-D-brane) system. Such a single-field effective action, to the lowest order in V^V^ around 
the top of the tachyon potential, can be obtained by the stringy computations for either a D3 
brane in bosonic theory |218ll219j or a non-BPS D3 brane in supersymmetric theory [220J (see also 
[22 lj for an application in inflation). What distinguishes the tachyon action from the standard 
Klein-Gordon form for a scalar field, is that the former is non-standard but it is of the "Dirac- 
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Born-Infeld" form [222.223.224j. The tachyon potential is derived from string theory itself and 
thus it has to satisfy some definite properties and requirements, such as to describe the tachyon 
condensation. 



6.2.1 A model inspired by string theory 

In this paragraph we construct a specific quintom model inspired by string theory. One considers 
an action of the form [110J: 



S = J d*Xy/=g~ [-V (</>)y/l - a'V M 0V^ + p'<f>n<f>\ , (6.22) 

which generalizes the usual "Born-Infeld-type" action, for the effective description of tachyon 
dynamics, by adding a term <pD<p to the usual V M 0V M in the square root. The two parameters 
a' and /3' in (I6.22p can be also made arbitrary when the background flux is turned on |225j . Note 
that we have defined a' = a/M 4 and /3' = /3/M 4 with a and /3 being the dimensionless parameters 
respectively, and M an energy scale used to make the "kinetic energy terms" dimensionless. 
V{4>) is the potential of scalar field <p (e.g., a tachyon) with dimension of [mass] 4 , with a usual 
behavior in general, that is bounded and reaching its minimum asymptotically. Note that, □ = 
~jTgdtiy/~99 >tv dvi an d therefore in ( 16.221) the terms V^0V M and </>□</> both involve two fields 
and two derivatives. 

The equation of motion for the scalar field <fi reads: 

§o(^) +aV „(rp) + ^ V+ f D ^0, ,6.23) 



where / = Jl — a'V^V^ + /?'</>□</> and V$ = dV/dcf). Following the convention of [87] , the 

energy-momentum tensor T^ v is given by the standard definition 8 g S = — J d i x 2 ^-T in '6g liV , 
thus: 

= gTWWf + V P (VV^)] + ^j^V>V> - V^V> - VW, (6.24) 
where if> = J§ = 

For a flat FRW universe and a homogenous scalar field <f>, the equation of motion (16.231) can be 
solved equivalently by the following two equations 

^ + 3 ^=^-f + ^, (6.25) 
V V 4M> 2 0(f) fi<y ' v ; 

4, + Mi, = (2a + „) - _) - -(2c - fl^ - ^ - , (6.26) 

where we have made use of if) as defined above. Therefore, (16.251) is just the defining equation 
for if) in terms of <f> and its derivatives. Finally, from (16. 24ft we obtain the energy density and 
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Fig. 17. (Color online) The DE EoS parameter as a function of Ina for V(<p) = V$e A< ^ 2 with Vq = 0.8, 
A = 1, a = 1, and j3 = —0.8, and for initial conditions <pi = 0.9, <pi = 0.6, (□(/>)« = ^( n <A)j = 0. From 
Ref. H22F. 

pressure: 



p = Vf + 



d 



P- 



-Vf- 



a 3 dt 
d 



o?dt 



(6.27) 
(6.28) 



The numerical exploration of the constructed model has been performed in [llOj . In Fig. [T7] 
one can see the EoS parameter of DE for the (motivated by string theory) potential V(<f>) = 
Voe'^ . The model experiences the phantom-divide crossing, and results to a Big Rip singularity. 
Numerically one sees that = when w crosses —1. 

In Figjl8]we assume a different potential, namely V(<p) = ^x$t^=x$ ■ One observes again that the 
EoS of the model crosses —1 during the evolution. Finally, in FigJT9]the potential has been also 
chosen as V(<f>) = pX4 ,+ P -\<t, ; but taking a positive j3 we observe that the EoS parameter starts 



e A0_j_ e — \<f> J 

with w < —1, it crosses —1 from below to above and then transits again to w < 



-1. 



One can generalize the aforementioned model in the case of a tachyon, non-minimally coupled 
to gravity [226] . Thus, one adds an extra term TOT to the usual terms in the square root and 
the non-minimal coupling Rf(T) |1191l 227j: 



S 



d x\ 



Ml 



?-Rf{T) - AV(T)Jl - a'g^d^TdJT + fi'TUT 



(6.29) 
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Fig. 18. (Color online) The DE EoS parameter as a function of Ina, for V(fa) = £ a^>^°-a0 with Vq = 0.5, 
X = 1, a = 1, and {3 = —0.8, and for initial conditions fa = 0.9, fa = 0.6, (d</>)i = ^( a< A)j = 0. From 

Ref. gnu. 
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Fig. 19. (Color online) The DE EoS parameter as a function of Ina, for V{fa) = e \^ e -\^ with Vq = 0.5, 
A = 1, a = 1, and /3 = 0.8, and for initial conditions fa = 0.9, fa = 0.6, (^fa)i = ^{^fyi = 0. From 

Ref msj. 

where M p is the reduced Planck mass. Action (I6.29f) can be brought to a simpler form for the 
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derivation of the equations of motion by performing a conformal transformation — > f(T)g llv : 



d x\/—g 



Ml 



3/ 



12 



2/ 



-^(R-^-d^Td fl T)-AV(T) x /l - (a'f(T) — 2/3'/ / (T)T)<9 A( T(9^T + /3'f(T)TaT 



(6.30) 

where V(T) = fp- is the effective potential of the tachyon. For a flat FRW universe and a 
homogenous scalar field T, the equations of motion can be solved equivalently by the following 
two equations, 



T + 3HT ■ 



'2/3' f'T - a'f_\^ f _ A 2 P'fVT ~„ 



fT 



2ip 



3M!(fff»-f* 



2 V f- 



rpZ 



(i -£')(«' -2/3') r 



a' 



f T 



y>r 2 



fTjfryf 2 - 2 la' - 2/3'^Tj Hf 



1 + f - 3f T 



3MI f\ T 
2 I / ) i> 



7 



(6.31) 



(6.32) 



where 



dC 
BUT 



A/3'VfT 
2h 



(6.33) 



h = y/l - (a'f - 2/3'f'T)d ll Td»T + pfTUT, and V 
pressure read [119IT2 27J: 



Finally, the energy density and 



p = AVh + 



d 



a 3 dt 



a 3 i)T ) + [a'f- 2(3'f'T ) ^-f 2 - 2^T + 



3M p 2 /f 



h 



7 17 



p = -AVh- 



d 

Jdt 



a 3 ^T + 



3M 2 ^ 
/ 



T 



(6.34) 
(6.35) 



The numerical exploration of the constructed model has been performed in |119j . In Fig. [201 we 
present the EoS parameter of DE for the (motivated by string theory) potential V(T) = e~ XT2 
and f(T) = 1 + Hi=\c{T 2% . Thus, the model experiences the — 1-crossing, before the tachyon 
potential reaches asymptotically to its minimum. 



6.2.2 Analysis on perturbations 

Let us examine the stability of the aforementioned model by considering quadratic perturbations. 
Consider a small perturbation 7r around the background (f>, 
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Fig. 20. (Color online) The plot of EoS parameter for the potential V(T) = e xt2 , for a = —2 and 
(3 = 2.2. The initial values are = 1 and = 3. From Ref. 1119]. 



0^0(t) + vr(£,x) , 



(6.36) 



where the background field in the FRW cosmology is spatially homogenous. After this shift of 
the field, and a tedious calculation, we obtain the following terms for the quadratic perturbations 
of the action |llUj : 



2 Vf + \ V1 ' 



-3 



7T 



(6.37) 



Interestingly, we notice that due to the positivity of the term ^-V / _3 2 in this model if the 
coefficient of (V7r) 2 is positive, the term in front of 7r 2 is guaranteed to be positively valued. The 
sound speed characterizing the stability property of the perturbations is given by 



a 
2 



vr 1 + f vf- 1 + W- 



-4>^4>vf- 



% Vf -i + f V f-^ + f Vf-' + f^/-i - ^0D0V/- 3 



(6.38) 



If the c 2 lies in the range between and 1 and the coefficient of 7r 2 remains positive, then the 
model will be stable. In Fig. [2TJ and Fig. [22] we present c 2 and the coefficients of 7r 2 respectively, 
for the models of Figs. [EJ dU and [HO We observe that for these models c 2 lies in the range [0, 1] 
and the coefficients of 7r 2 are positive, thus the models are indeed stable. 
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Fig. 21. (Color online) Evolution of the sound speeds (from left to right) for the three models considered 
in the text. From Ref. lllOf . 
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Fig. 22. (Color online) Evolution of the coefficient of-k 2 (from left to right) for the three models considered 
in the text. From Ref. Ill Of . 
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7 Realizations of quintom scenario in non-scalar systems 



7.1 Spinor quintom 



The quintom models described so far have been constructed by scalar fields, which are able 
to accommodate a rich variety of phenomenological behaviors, but the ghost field may lead 
to quantum instabilities. To solve this problem one may consider the linearized perturbations 
in the quantum-corrected effective field equation at two- loop order [228(2 29, 230J and obtain 
a super-acceleration phase induced by these quantum effects [231 |2 32j. However, there is the 
alternative possibility where the acceleration of the universe is driven by a classical homogeneous 
spinor field. Some earlier studies on applications of spinor fields in cosmology are given in Refs. 
|233|[2 34 235j. In recent years there are many works on studying spinor fields as gravitational 
sources in cosmology. For example see Refs. [236|237|238j for inflation and cyclic universe driven 
by spinor fields, Refs. [239] for spinor matter in Bianchi Type I spacetime, Refs. |240|241j for a 
DE model with spinor matter etc. 



7.1.1 Algebra of a spinor field in FRW cosmology 

Recently, a delicate quintom model was proposed by making use of spinor matter [122] . This 
scenario can realize the quintom behavior, with an EoS crossing —1, without introducing a ghost 
field. Instead, the derivative of its potential with respect to the scalar bilinear tpip, which is 
defined as the mass term, becomes negative when the spinor field lies in the phantom-like phase. 
If this model can realize its EoS across —1 more than one time, the total EoS of the universe can 
satisfy w > — 1 during the whole evolution, as it is required by NEC [242f243] . and we expect to 
treat this process as a phase transition merely existing for a short while. 

To begin with, we simply review the dynamics of a spinor field which is minimally coupled 
to Einstein gravity (see Refs. [244,245,246] for a detailed introduction). Following the general 
covariance principle, a connection between the metric g^ and the vierbein is given by 

9^a e l = Vab , (7.1) 

where denotes the vierbein, g^ is the space-time metric, and r] a b is the Minkowski metric with 
7] a t, = diag(l, —1, —1, —1). Note that the Latin indices represents the local inertial frame, while 
the Greek indices represents the space-time frame. 

We choose the Dirac-Pauli representation as 



(i o\ 




( a\ 






K o-i) 







(7.2) 



where <Ji is Pauli matrices. One can see that the 4x4 7° satisfy the Clifford algebra{7 a , 7 fe } = 2r\ a ^. 
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The 7° and provide the definition of a new set of Gamma matrices 



F" = e ^ a , (7.3) 

which satisfy the algebra T^} = 2g fJLU . The generators of the Spinor representation of the 
•entz group can be written as S ab = 

by 



Lorentz group can be written as S ab = i[7 a ,7 b ]. Therefore, the covariant derivatives are given 



(7.4) 
(7.5) 



where the Dirac adjoint ip is defined as ip + ^° . The 4x4 matrix f2 M = ^u^ a bT l ab is the spin 
connection, where oj^ab = e^V^e^ are the Ricci spin coefficients. 

Using the aforementioned algebra we can write the following Dirac action in a curved space-time 
background 



d x e 



(7.6) 



Here, e is the determinant of the vierbein and V stands for the potential of the spinor field 
if) and its adjoint ip. Due to the covariance requirement, the potential V depends only on the 
scalar bilinear ipip and on the "pseudo-scalar" term tpj^ip. For simplicity we neglect the latter 
term and assume that V = V(if)if)). 

Varying the action with respect to the vierbein e^, we obtain the energy-momentum-tensor 



T, 



e 5e u „ 



(7.7) 



On the other hand, variation of the action with respect to the fields ip, ip respectively yields the 
equation of motion of the spinor 



dV - dV 



(7- 



As usual, we deal with the homogeneous and isotropic FRW metric. Correspondingly, the vierbein 
are given by 



a 



(7.9) 
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Assuming that the spinor field is space-independent, the equations of motion read: 



ij) + -HiJ) + i'fV'il> = (7.10) 
i + -H$ -i 7 VV = , (7.11) 

where a dot denotes a time derivative £ 4', a prime denotes a derivative with respect to t/^, and 
if is the Hubble parameter. Taking a further derivative, we can obtain: 



N 

•W> = -3 , (7-12) 
where iV is a positive time-independent constant, defined as the value of ipip at present. 

From the expression of the energy- momentum tensor in (17. 7ft . we obtain the energy density and 
pressure of the spinor field: 



Pi , = T° = V (7.13) 
P4, = -T! = V'W>-V , (7.14) 

where Eqs. (I7.10p and (17. lip have been applied. The EoS of the spinor field, defined as the ratio 
of its pressure to energy density, is given by 



w lf) = — = -l-\ r— . (7.15) 



If we assume the potential to be a power law as V = V (^) a , with a as a nonzero constant, we 
obtain a constant EoS: 



itty, = — 1 + a . (7-16) 

In this case, the spinor matter behaves like a linear-barotropic-like perfect fluid. For example, 
if a = ~, we can acquire ~ a -4 and uy, = ~, which is the same as radiation. If a — 1, then 
~ a~ 3 and = 0, thus this component behaves like normal matter. 

It is interesting to see that the spinor matter is able to realize the acceleration of the universe 
if a < |. Thus, it provides a potential motivation to construct a dynamical DE model with the 
spinor matter. In the next subsection we focus our investigation on constructing the subclass of 
quintom DE models using the spinor field, which is called Spinor Quintom. 
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7.1.2 Evolutions of Spinor Quintom 

In the model of the previous subsection V is positive in order for the energy density to be positive 
(see (I7.13P ). Therefore, from (17.151) we deduce that > —1 when V > and < —1 when 
V < 0. The former corresponds to a quintessence-like phase and the latter to a phantom-like 
phase. Therefore, the potential derivative V is required to change its sign if one desires to obtain 
a quintom picture. In terms of the variation of V there are three categories of evolutions in 
spinor quintom: 

(i) 



V > o -> V < o , 

which gives a Quintom-A scenario by describing the universe evolving from quintessence-like 
phase with > — 1 to phantom-like phase with < — 1; 



V < o -> 7' > o , 

which gives a Quintom-B scenario for which the EoS is arranged to change from below —1 to 
above — 1; 

(iii), V changes its sign more than one time, thus one obtains a new Quintom scenario with its 
EoS crossing —1 more than one time, dubbed Quintom-C scenario. 

In the following, we assume different potentials in order to realize the three types of evolutions 
mentioned above GO- 

To begin with, we shall investigate Case (i) and provide a Quintom-A model. We use the potential 
V = Vq[{i[)iP — 6) 2 + c], where V , b, c are undefined parameters. Then we obtain V = 2V (ijjip — b) 
and the EoS reads: 



~ fe 2 ~ C 

* {^) 2 - 2hfjtfj + b 2 + c' 1 ' ' 

According to (17. 12ft one finds that ipip is decreasing with increasing scale factor a during the 
expansion of the universe. From the formula of V we deduce that at the beginning of the 
evolution the scale factor a is very small, so t/np becomes very large and ensures that initially 
V > 0. Then tpip decreases with expanding a. At the moment where = b, one can see 
that V = 0, which results in the EoS — — 1. After that V becomes less than and thus the 
universe enters a phantom-like phase. Finally, the universe approaches a de-Sitter spacetime. This 
behavior is also obtained by numerical elaboration and it is shown in Fig. [23 From this figure 
one observes that the EoS w^p starts the evolution from 1, then mildly increases to a maximum 
and then begins to decrease. When ipip = b, it reaches the point w^p = — 1 and crosses —1 from 
above to below smoothly. After that, the EoS progressively decreases to a minimal value, then 



3 Note that we choose the potentials phenomenologically without any constraints from quantum field 
theory or other consensus. From the phenomenological viewpoint it is allowed to treat the background 
classically, while dealing with the perturbations in quantum level, like what it is usual in inflation theory. 
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Fig. 23. (Color online) Evolution of the EoS in Case (i) as a function of time. For numerical elab- 
oration we consider the potential of the spinor field as V = VoK^ip — b) 2 + c], where we choose 



Vq = 1.0909 x 10 , b = 0.05 and c = 10 for the model parameters. For the initial condition 
we take N = 0.051. From Ref. fW% . 

increases and eventually approaches the cosmological constant boundary. 



In Case (ii) we choose the potential as V 
Vo(-2$V + b ) and the EoS 



Vq[— {ipip — b)ipip + c]. Then one can obtain V 



— (vjjijj) 2 + bipif) + c 



(7.18) 



Initially V is negative due to the large values of ipip. Then it increases to when ipip = |, 
whereafter it changes its sign and becomes larger than 0, in correspondence with the Case (ii). 
From Fig. [2U we can see that the EoS evolves from below —1 to above —1, and finally it 
approaches —1. 

In the third case we explore a quintom scenario which presents the — 1-crossing for two times, 
using the potential V = Vo[(ijjip ~ b) 2 vbip + c]. Thus, we acquire V = V$(ipip — b)(3tjjip — b) and 
the EoS becomes 



2(^) 3 - 26(^) 2 - c 
(^) 3 - WfW + b 2 {i)yb) + c 



(7.19) 



Evidently, the equation V = has two solutions which are = b and = |, thus V changes 
its sign two times. From the EoS expression (17.191) . we find that initially uy, > —1. When the 
value of tpifj becomes equal to b, it crosses —1 for the first time. After the first crossing, it 
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Fig. 24. (Color online) Evolution of the EoS in Case (ii) as a function of time. For numerical elab- 
oration we consider the potential of the spinor field as V = Vo[— (ifiip — b)ij)tp + c], where we choose 
Vo = 1.0909 x 10~ 117 , b = 0.05 and c = 10~ 3 for the model parameters. For the initial condition we take 
N = 0.051. From Ref. f!2^ . 

enters the phantom-like state, it continuously descends until it passes through its minimum, 
then it ascends to Tpip = | and then it experiences the second crossing, moving eventually to the 
quintessence-like phase. This behavior is depicted in Fig. One can see that in this case the 
Big Rip can be avoided. 

Moreover, considering the dark-matter component with energy density pu oc 1/a 3 and % = 0, 
we can see that the EoS of the universe w u = p P + PM satisfies the relation w u > — 1 in this 
case. As it is argued in Ref. |243] . NEC might be satisfied in such a model, since < —1 is 
only obtained for a short period during the evolution of the universe. 

As mentioned in previous sections, one important issue of a DE model is the analysis of its 
perturbations. Performing such an investigation we might be able to learn to what degree the 
system is stable both at quantum and classical level. Usually systems with w < — 1 present 
instabilities, for example a phantom universe suffers from a Big Rip singularity. Although this 
classical singularity can be avoided in most quintom models, which usually enter a de-Sitter 
expansion in the final epoch, all the scalar quintom models suffer from quantum instabilities 
since there are negative kinetic modes from ghost fields. Here we would like to perform the 
perturbation theory of spinor quintom. Since we do not introduce any ghost fields in our model, 
the phantom divide crossing is achieved by the spinor field itself and thus it does not result to 
any particular instability. 

3 

In order to simplify the derivation, we would like to redefine the spinor as ip^ = azip. Then 
perturbing the spinor field we extract the perturbation equation as 
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Fig. 25. (Color online) Evolution of the EoS in Case (Hi) as a function of time. The magenta solid line 
stands for the EoS of spinor quintom and the green dot line stands for that of the whole universe. For 
numerical elaboration we consider the potential of the spinor field as V = VoKtptp — b) 2 ^ 1 ^ + c], where we 
choose Vq = 1.0909 x 10~ 117 ; b = 0.05 and c = 10~ 3 for the model parameters. For the initial condition 
we take N = 0.051. From Ref. fW% . 
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(7.20) 



where r is the conformal time defined by dr = dt/a. Since the right hand side of the equation 
decays proportionally to a~ 3 or even faster, we can neglect those terms during the late-time 
evolution of the universe for simplicity. 



From the perturbation equation above, we can read that the sound speed is equal to 1, and this 
eliminates the instability of the system in short wavelength. Furthermore, when the EoS w crosses 
— 1 we have V = and thus the eigenfunction of the solution to Eq. (I7.20p in momentum space 
is a Hankel function with an index |. Therefore, the perturbations of the spinor field oscillate 
inside the hubble radius. This is an interesting result, because in this way we might be able to 
establish the quantum theory of the spinor perturbations, just as what it is done in inflation 
theory. 



Note that the above derivation does not mean that spinor quintom is able to avoid all instabilities. 
We still have not examined the effects of the right hand side of Eq. (I7.20p . which could destroy 
the system stability under particular occasions. Finally, another possible instability could arise 
from quantum effects, if this model is non-renormalizable. 
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7.1.3 A unified model of Quintom and Chaplygin gas 

As we have analyzed so far, a spinor field with a power-law-like potential behaves like a perfect 
fluid with a constant EoS. However, it is still obscure to establish a concrete scenario to explain 
how a universe dominated by matter evolves to the current stage, that is dominated by DE. 
The last years, an alternative interesting perfect fluid, with an exotic EoS p = —c/p, has been 
applied into cosmology |248] in the aim of unifying a matter-dominated phase where p oc 1/a 3 , 
and a de-Sitter phase where p = —p, which describes the transition from a universe filled with 
dust-like matter to an exponentially expanding universe. This so-called Chaplygin gas [ 248] and 
its generalization [249] have been intensively studied in the literature. 

Some possibilities of this scenario, motivated by field theory, are investigated in [250] . Inter- 
estingly, a Chaplygin gas model can be viewed as an effective fluid associated with D-branes 
|251|l252 j. and it can be also obtained from the Born-Infeld action |253|254j . The combination of 
quintom and Chaplygin gas has been realized by the interacting Chaplygin gas scenario [255], 
as well as in the framework of Randall-Sundrum braneworld [256] . Finally, the Chaplygin gas 
model has been thoroughly investigated for its impact on the cosmic expansion history. A con- 
siderable range of models was found to be consistent with SNIa data [257] , the CMBR [258] , the 
gamma-ray bursts [259], the X-ray gas mass fraction of clusters [26U] . the large-scale structure 
[2BTj . etc. 

Here, we propose a new model constructed by spinor quintom, which combines the property of 
a Chaplygin gas. The generic expression of the potential is given by 

V = + c , (7.21) 

where f{$ip) is an arbitrary function of Altering the form of f(ipifj), one can realize both 
the Chaplygin gas and quintom scenario in a spinor field. 

Firstly, let us show how this model recovers a picture of generalized Chaplygin gas. Assuming 
f(ipip) = Vq (?/>■?/>) 1+/3 the potential becomes 



V = ^o(#)^ + c , (7.22) 
and therefore, the energy density and pressure of (17. 131) . (17. 141) read 

P^= 1+ {/v (^y+? + c (7.23) 

Pi, = -c[V (^) 1+/? + c]~W . (7.24) 



Thus, this spinor model behaves like a generalized Chaplygin fluid, satisfying the exotic relation 
Pip — —~r- Being more explicit we assume /3 = 1, obtaining the expressions of energy density 

and pressure as 
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^ = \/-6- + C > Pl> = -— ■ ( 7 - 25 ) 



In this case a perfect fluid of Chaplygin gas is given by a spinor field. Based on the above analysis, 
we conclude that this simple and elegant model is able to mimic different behaviors of a perfect 
fluid and therefore it accommodates a large variety of evolutions phenomenologically. 

In succession, we will use this model to realize a combination of a Chaplygin gas and a Quintom- 
A model, which is mildly favored by observations. Choosing fi^ip) to be f(?p?p) = Vq(^i/; — b) 2 
we acquire the potential 



V = vV (V^ -b) 2 + c , (7.26) 
where Vq, b, c are undetermined parameters. Its derivative reads 

V (M-b) 

V = '= , (7.27) 

yVo^ -bf + c 

and thus the EoS becomes 

with the — 1-crossing taking place when ipijj = b. During the expansion ipifj is decreasing following 
f)7.12p . with ipifj — > oo at the beginning of the evolution where a — > 0. Therefore, increases 
from to a maximum value and then it starts to decrease. When {/jip = b the EoS arrives at the 
cosmo logical constant boundary uy, = — 1 and then it crosses it. Due to the existence of c-term, 
the EoS eventually approaches the cosmological constant boundary. In this case the universe 
finally becomes a de-Sitter space-time. 

Considering a universe filled with dark matter and the aforementioned spinor quintom matter, we 
perform a numerical analysis and we present the evolution of the EoS in Fig. [2(3 Moreover, in Fig. 
[27] we display the evolution of the density parameters of dark matter Q M = Pm/(pm + P^) and 
spinor quintom D,^ = p^j (pm +Pip)- It is evident that the ratio of these two components from the 
beginning of evolution is of order one, in relief of the fine-tuning and coincidence problems. During 
the evolution DE overtakes dark matter, driving the universe into an accelerating expansion at 
present, and eventually it dominates the universe completely, leading to an asymptotic de-Sitter 
expansion. 

Summarizing the results of this section we see that while a scalar-type quintom model leads to 
quantum instabilities (due to the ghost field with a negative kinetic term), the class of quintom 
models in virtue of a spinor field can avoid such a problem. As shown above, such a scenario could 
behave as a generalized Chaplygin gas, and give rise to —1-crossing during the evolution, caused 
by the sign-flip of the V'-term. Compared with other models experiencing the —1-crossing, this 
scenario is also economical, in the sense that it involves a single spinor field. 
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Fig. 26. (Color online) Evolution of the EoS of the unified model with potential \7.26 ) as a function of 
time. For the numerical analysis we assume Vq = 3.0909 x 10~ 239 , b = 0.05 and c = 9 x 10 -241 , while 
for the initial conditions we take N = 0.051. From Ref. 1122]. 



1 1 1 1 1 1 1 1 1 1 1 1 1 










Q 












Q 






m 




\ / 






\ / 






/ \ 

■ — y \ 






\ 






\ 






\ 






\ 












\ 

\ 






\ 

\ 






1 . 1 . 1 . 1 . 1 


i 


i 



-4 -2 2 4 6 8 



t 



Fig. 27. (Color online) Evolution of the density parameters of DE (violet solid line) and dark matter 
(red dash line) as a function of time. From Ref. [122]. 
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7.2 Other non-scalar models 

We finish this section mentioning that non-scalar systems may give rise to a quintom scenario 
under certain cases. For example, developed from a vector-like DE model [262 263J, Ref. [264J 
has shown that an unconventional cosmological vector field in absence of gauge symmetry can 
realize its EoS across —1. This model and its extensions have also been shown to be possi- 
ble candidates for DE |265|I266|I267||2"6 8 269 2 70f271j . and even solve the coincidence problem 
|272f273f274j . However, see |275] for a different viewpoint. 

Finally, recently it has been noticed that the quintom scenario can be achieved in a model of non- 
relativistic gravity [276|I277|I278J . but this is realized at a high energy scale, far away from that 
of DE. The possible successes of this model has led to the exploration of the potential relations 
among quintom scenario and phenomena of Lorentz symmetry breaking [279,280,281,282,283]. 




8 Quintom scenario in the braneworld 

An alternative way of explaining the observed acceleration of the late universe is to modify grav- 
ity at large scales. A well-studied model of modified gravity is the braneworld model. Although 
the exciting idea that we live in a fundamentally higher-dimensional spacetime which is greatly 
curved by vacuum energy was older [28412850286028702880289029002910292029302940295] . the new 
class of "warped" geometries offered a simple way of localizing the low energy gravitons on the 
4D submanifold (brane) 296.2970298) . Motivated by string/M theory, AdS/CFT correspondence, 
and the hierarchy problem of particle physics, braneworld models were studied actively in re- 
cent years [299,300,301,302,303,304,305], and have been shown to admit a much wider range of 




possibilities for dark energy |306j . A well studied model is the Dvali-Gabadadze-Porrati (DGP) 
one [307,308j where the braneworld is embedded in the flat bulk with infinite extra dimension. 
In this model, gravity appears 4D at short distances but it is altered at distance larger than 
some freely adjustable crossover scale r c , through the slow evaporation of the graviton off our 
4D braneworld universe into the fifth dimension. The inclusion of a graviton kinetic term on the 
brane recovers the usual gravitational force law scaling, 1/r 2 , at short distances, but at large 
distances it asymptotes to the 5D scaling, 1/r 3 . The matter particles cannot freely propagate in 
the extra dimensions, and are constrained to live on the (tensionless) brane. In such a model, 
late-time self-acceleration solutions appears naturally [309f310Q311j . since they are driven by the 
manifestation of the excruciatingly slow leakage of gravity into the extra dimension. 

8.1 Quintom DE in DGP scenario 

In this subsection we investigate a simple model of a single scalar field in the framework of DGP 
braneworld, following [108] . Although the behavior of the effective DE on a DGP brane with a 
cosmological constant and dust has been studied in [10203 12] . from the unified theoretic point of 
view we can assume that the gravitational action is not necessarily the Einstein-Hilbert action. 
In fact, string theory suggests that the dimensionally reduced effective action includes not only 
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higher-order curvature terms but also dilatonic gravitational scalar fields. Thus, at the level of 
the low-energy 5D theory (we restrict to 5D braneworld models although even higher- dimensional 
are also possible |313f314f315f31(jf317] ). it is naturally expected that there appears a dilaton-like 
scalar field in addition to the Einstein-Hilbert action |318. 319f32U] . Hence, it is interesting to 
investigate how such a scalar field in the 5D theory affects the braneworld |321|322j . Since we 
consider a single-field model, the accelerated expansion of the universe will be a result of the 
combined effect of the field evolution and the competition between 4D gravity and 5D gravity. 

Let us start from the action of the DGP model 



S — SVmlk + ^ 



'brane; 



where 



and 



S 



bulk 



(FX 



M 



-(5), 



2k§ 



•Sbrane — Id X\J ~Q 
JM 



K 5 



+ L 



branc 



$.2) 



•3) 



Here k\ is the 5-dimensional gravitational constant, is the 5-dimensional curvature scalar, 
x M (/i = 0, 1, 2, 3) are the induced 4D coordinates on the brane and K is the trace of extrinsic 
curvature on either side of the brane. Finally, L hra _ ne (g a /3, if)) is the effective 4D Lagrangian, which 
is given by a generic functional of the brane metric g a p and brane matter fields if>. 



We consider a brane Lagrangian consisting of the following terms 

,2 



P 



7 brane 



R + L m + L a 



•4) 



where p is the 4D reduced Planck mass, R denotes the curvature scalar on the brane, L m stands 
for the Lagrangian of other matters on the brane, and represents the lagrangian of a scalar 
confined to the brane. Assuming a mirror symmetry in the bulk we obtain the Friedmann equation 
on the brane 



H 2 + 



1 

3pJ 



p + Po + 6po(l + 



2p\V 2 



Po 



i.5) 



where k is the spatial curvature of the three dimensional maximally symmetric space in the FRW 
metric on the brane, and 9 — ±1 denotes the two branches of DGP model, p denotes the total 
energy density on the brane, including dark matter and the scalar field (p = p^ + pdm), and the 
term p relates the strength of the 5D gravity with respect to 4D gravity, that is p = where 
the cross radius is defined Comparing the modified Friedmann equation on the 



brane with the standard one: 



H 2 + 



k 



Pdm + Pdt 



1 

a* 3/i 2 

one obtains the density of the effective 4D dark energy as 



Po + 0p c 



p + po + 0pa[l + 



2 P y/ 2 

Po 



1.6) 



i.7) 
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As usual DE satisfies the continuity equation: 

dpde 



dt 



+ 3H(p de +p e ff) = 0, 



where p e ff denotes the effective pressure of DE. Then we can express the EoS of DE as 

Peff ,1 d\np de 
Wde = = -1 + -- 



Pde 



3dln(l + z) 



.9) 



where, from (18. 8 \ 

d In p de , 3 



dln(l + z) pde 



Pde +Pde + 9 



1 + 2 



Pde + Pdr, 



-1/2 



+ Pdm + Pde 



.10) 



We mention that the DE pressure p de is different from p. 



eff- 



Clearly, if p de decreases and then increases with respect to redshift, or increases and then de- 
creases, it is implied that DE experiences a phantom-divide crossing. Equations like (I8.5P and 
(18.81) coincide with the subclass of inhomogeneous EoS of FRW universe [90] or FRW universe 
with general EoS [323J, which can give rise to the phantom-divide crossing. 



8.1.1 Canonical scalar field 



For a canonical scalar with an exponential potential of the form: V = V^e '< the effective DE 
evolves as [108J: 



dpde 



rfln(l + z) 



+ 1 + 



+ 2V + 2p 



dm 



-1/2 



+ Pdr, 



(8.11) 



If 8 = 1 both terms on the RHS are positive and thus it never goes to zero at finite time, while in 
the case 9 = — 1 the two terms of the RHS have opposite sign, therefore it is possible to obtain 
the — 1-crossing. Hence, we restrict our study to 9 = — 1. 

In order to transform the aforementioned dynamical system into its autonomous form 
we introduce the auxiliary variables: 



x 



I 



>Pm 



\^6pH ' \f3pH ' \/3pH ' ^/3pH 

Thus, the dynamics of the system can be described by: 



(8.12) 
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x' = ~ax{2x 2 + l 2 ) + 3x-^j-Xy 2 , (8.13) 

y ' = J-ay{2x 2 + l 2 ) + ^-\xy 1 (8.14) 
l' = --al(2x 2 + / 2 ) + -/, (8.15) 



2 v ' 2 

3 
— ( 

2 



6' = -^6(2x 2 + / 2 ), (8.16) 



where a = 1 — (1 + 2 X +l )~ 1//2 , and here a prime stands for derivation with respect to 
s = — ln(l + ,z), and we have set k — 0. One can see that this system degenerates to a quintessence 
with dust matter in standard general relativity. Note that the four equations (18.131) . (18.141) . (18.151) . 
(18.161) are not independent, since the constraint 



„2 i „,2 i ;2\ V 2 



x * + y* + P + b 2 -b 2 h+2 X +V b2 +l 1 =1, (8.17) 

arising from the Friedmann equation, reduces the number of independent equations to three. 
There are two critical points of this system, satisfying x' = y' = I' = b' = 0, appearing at 



x = y = I = 0, b = constant; (8.18) 
x=y = 1 = b = 0. (8.19) 

However, neither of them satisfies the Friedmann constraint (18.171) . Hence one proves that there 
is no (kinetic energy)- (potential energy) scaling solution or (kinetic energy)- (potential energy)- 
(dust matter) scaling solution on a DGP brane with quintessence and dust matter. 

The most significant parameter from the viewpoint of observations is the deceleration parameter 
q, which carries the total effects of cosmic fluids and it is defined as q = —™ = — 1 + ^a(2x 2 + l 2 ). 
In addition, for convenience we introduce the dimensionless density as 



p ^Pde = SK f 2 + y2 + b 2_ b 2 

Or b 2 



1 + 2 



x 2 + y 2 + I' 
b 2 



1/2' 



(8.20) 



where p c denotes the present critical density of the universe, and the rate of change with respect 
to redshift of DE: 



In Fig. [28] we present a concrete numerical example of the — 1-crossing, depicting 0, 7 and q as 
a function of s = — ln(l + z). Thus, contrary to the conventional 4D cosmology, where a single 
canonical field cannot experience the —1-crossing, in the present model this is possible due to 
the induced term p of the "energy density" of r c . Only a small component of po> i- e - ^r c — 0.01, 
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Fig. 28. (a) /3 (solid curve) and 7 (dotted curve) as functions of s = — ln(l + z) for the canonical field 
case of the DGP model. The EoS parameter of DE crosses — 1 at about s = —1.6, i.e. at z = 3.9. (b) 
The corresponding deceleration parameter q vs s, which crosses at about s = —0.52, i.e. at z = 0.68. 
For this figure, VL^ = 0.01, VLm = 0.3, and VL Tc = po/p c = 0.01, A = 0.05. From Ref. \108l. 

is capable of making the EoS of DE cross —1. At the same time the deceleration parameter is 
consistent with observations. 



8.1.2 Phantom scalar field 

For a phantom scalar with an exponential potential of the form: V = V$e 1* the effective DE 
evolves as [108J: 



dp, 



de 



dln(l + z) 



1 



-<P 2 + 2V + 2p, 



dm 



-1/2 



-<P + Pdn 



(8.22) 



If 9 = —1, both terms of RHS are negative and it never goes to zero at finite time. Contrarily, 
if 9 = 1 the two terms of RHS have opposite sign and the phantom-divide crossing is possible. 
In the following we consider the branch of 9 = 1. The dynamics is described by the following 
autonomous system: 



x' = -^ a 'x(-2x 2 + I 2 ) + 3x- ^-Xy 2 , (8.23) 

y> = -\a'y{-2x 2 + I 2 ) + ^Xxy, (8.24) 
l' = -^ a 'l(-2x 2 + / 2 ) + i, (8.25) 
b' = -^ a 'b(-2x 2 + / 2 ), (8.26) 

where a' = 1 + (1 + 2 ~ x ffi +l )~ 1/ ' 2 , and we have also adopted the spatial flatness condition. Now 
the Friedmann constraint becomes 
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Fig. 29. (a) /3 (solid curve) and 7 (dotted curve) as functions of s = — ln(l + z) for the phantom field 
case of the DGP model. The EoS parameter of DE crosses — 1 at about s = —1.25, i.e. at z = 1.49. (b) 
The corresponding deceleration parameter q vs s, which crosses at about s = —0.50, i.e. at z = 0.65. 
For this figure, Q ki = 0.01, Q M = 0.3, and n r - c = p /p c = 0.01, A = 0.01. From Ref. JIM) . 

- x 2 + y 2 + I 2 + b 2 + b 2 (l + 2 ~ X2+ J 2 + /2 ) V2 = 1, (8.27) 

with which there are three independent equations left in the system. Through a similar analysis as 
in the case of a canonical scalar field, it can be shown that there is no (kinetic energy)- (potential 
energy) scaling solution or (kinetic energy)- (potential energy)- (dust matter) scaling solution on 
a DGP brane with phantom and dust dark metter |108j . The deceleration parameter q becomes 
q = — 1 + \ct'(— 2x 2 + I 2 ), and the dimensionless density and rate of change with respect to 
redshift of DE become, 



P = ^l-x 2 + y 2 + b 2 + b 2 



1 + 2 



_ x 2 +2/ 2 + / 2, 1 l/2 



7 = 3 



-x 2 + y 2 + l 2 \]-V 2 
b 2 



I) 2 



'-2x 2 + I 2 ) + 2x 2 



(8.28) 
(8.29) 



In Fig. [29l we depict /3, 7 and q as a function of s = — ln(l + z). As we observe, the (effective) 
EoS parameter of DE crosses —1 as expected. This is also in contrast with conventional 4D 
cosmology, where a single phantom field lies always below —1. The 5D gravity plays a critical 
role in the realization of the — 1-crossing, and at the same time the deceleration parameter is 
consistent with observations. 



Finally, we stress that the most important difference of the dynamics between an ordinary and 
a phantom field is that while for the canonical field the —1-crossing takes place from above to 
below, in the phantom case it takes place from below to above (see 7 vs s in Figs. 1281 and 1291) . 



8.2 Quintom DE in the braneworld 



Let us now investigate the effects of the bulk quintom field in the DGP braneworld scenario, and 
in particular whether it is possible or not to have a late-time accelerated phase on the brane for 
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i) a tensionless brane with a quintom fluid in the bulk and ii) a quintom DE fluid on the brane 
and an empty bulk. We consider a 5D metric of the form 



-n 2 (t,y)dt 2 + a 2 (t,y)'jijdx 1 'dx j + b 2 (t,y)dy 2 



(8.30) 



where y is the coordinate of the fifth dimension and 7^ is a maximally symmetric 3-dimensional 
metric. We will use k to parameterize the spatial curvature and assume that the brane is a 
hypersurface located at y = 0. We shall be interested in the model described by the action 



S = J d 5 xV^g{J^R {5) - A + CT) + / d 4 x^(-5 b + C™ at ), 



(8.31) 



where is the curvature scalar of the 5D metric qab-, A is the bulk cosmo logical constant, 5 b 
is the brane tension, k\ = 8HG5, and qAB = Qab — n^ns {jia is the unit vector normal to the 
brane and A,B = 0,1,2,3,5) is the induced metric on the 3-brane. The 5D Einstein equation 
can be written as 



Rab - ^9abR {5) = i4[-Ag AB + Tab + S^5 B 5(y b )]. 



1 



(8.32) 



Here 5(y b ) = ^p-, Tab is the energy momentum tensor of the bulk matter, and the last term 
corresponds to the matter content in the brane 
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(8.34) 



(8.35) 
(8.36) 
(8.37) 



where primes indicate derivatives with respect to y and dots derivatives with respect to t. On 
the brane we assume a perfect fluid 



rj? = diag(-p b ,p b ,p b ,p b ). 



(8.38) 
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In the bulk space we consider a quintom field, constituted from the normal scalar field (p(t, y) 
and the negative-kinetic scalar field a(t,y), with a Lagrangian of the form: 



1 



pmat x AB 

~2 9 



; A <t>,B ~ v,a°,b J + V(<f>, a). (8.39) 
According to this action, the energy momentum tensor of the bulk quintom field is given by 



T 



AB 



;a<P,b - &,A a ,B - 9ab 



1 



.CD I 



,d ~ v,c<t,d) + V(4>,a) 



(8.40) 



Finally, the equations of motion of the scalar fields <f> and cr in the bulk space write: 



-a - 3-cr + ( n + 3— ) a' + a" + 

a V a / oa 



a 



n , dV{<f>,a) 



S£ mat 
Sa 



IA1) 
5.42) 



We are interested in studying the Einstein equations in the presence of a bulk quintom field at the 
location of the brane. Without loosing generality we choose b(t,y) = 1 and n(t, 0) = 1, which can 
be achieved by scaling the time coordinate. As is well known, the presence of the brane leads to a 
singular term proportional to 5-function in y on the right-hand sides of the Einstein equations (|8.34l) 
and (|8.35p and the equation of motions (|8.4ip . (j8.42p . which have to be matched by singularity in the 
second derivatives in y on the left-hand side. Since all fields under consideration are symmetric under 
the orbifold symmetry Z2, these jumps in the first derivatives in y fix the first derivatives at y = 0. In 
our case, the junction conditions read 



a' k| / 

— \y=0 = jr[Pb + Ob 

a 6 V 

n'\ y= o = y [3p b + 2p b - 5 b 



1A3) 



and 



\y=o 



1 bOl mt 



7 C \y=0 



1 bOl mt 



2 6<j) 1 ' Iy_v 2 Sa 
Using the 00 and 55 components of the Einstein equations in the bulk, one obtains 



.44) 



F' = ^(A-T °)aV-Mr 5 Vd, 
^ = ^(A-T|)a 3 d-M n2r o a 3^ 

where F is a function of t and y defined by 



(8.45) 
(8.46) 
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F (t,y) = ^>-( a 'a) 2 + ka 2 . (8.47) 

In order to find the solution for F, we assume that the quintom field in the bulk is independent of y. 
Under this assumption the non-vanishing components of the quintom energy momentum tensor take 
the form 



n = -p B = --^ + -a 2 -V{^a) 

n=p B = ^ 2 -\a 2 -V{<f>,a) 
T05 = W (re' = 

Tl= l -tf-U 2 -V{M, (8.48) 

valid also on the brane. As we observe, in the case of y-independent bulk quintom field, Tqs vanishes, 
thus there is no matter-flow along the fifth dimension. Furthermore, it is obvious from (|8.44|) that the 
quintom field cannot appear in the matter content on the brane |3211 322j. 

Solving equation (|8.45p leads to the first integral of the 00 component of Einstein equation as 



k\ { . \ C a? a' 2 k , 
^( A + M ) + ___ + ___ = , (8.49) 

where C is a constant of integration usually referred as dark radiation [337J. Substituting the junc- 
tion conditions f|8.43j) into (|8.49p . we arrive at the generalized Friedmann equation on the brane as 
|324l325|326f32 7 328 329l330|331|332|333l334l335|336|337|338| : 



* 2 4 - ! ( a + 1*) + 1** + i-> + P 4- « 8 -»» 

As one can see from (j8.50|) . in the absence of the bulk matter field, cosmological constant and tension, 
the equation gives rise to a Friedmann equation of the form H oc p\j instead of H oc yfp~b, which 
is inconsistent with cosmological observations. This problem can be solved by either considering the 
cosmological constant and tension on the brane or considering a matter field in the bulk |339|340j . 
Recalling the junction conditions (|8,43p . the 05 component of Einstein equations and the field equation 
(18.411) . (I8.42P acquire the following forms on the brane: 



p b + 3H(p b +p b ) = , (8.51) 
d\ . dV(4>, a 



a + 3(-)a ^^ = 0. (8.53) 



It should be noted that if the scalar fields <j> and a satisfy the field equations (|8.52p and (|8.53p respectively, 
the bulk energy momentum tensor is conserved and we obtain 
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Pb + SH( Pb +pb) = 0- 



(8.54) 



We are interested in studying the acceleration condition for a universe with quintom field in the bulk. 
The condition for acceleration can be obtained from (|8.5Up by using the conservation equation of the 
brane and bulk matter ( (|8.5ip and (|8.54p respectively): 



| = |(A + - g*(* + *) - |(,B + 3 PB ) - g(« + ta ) - £. (8.55) 

In order to study the role of bulk quintom field in the late-time acceleration phase on the brane, we 
neglect the effect of tension, brane matter and dark radiation, and therefore the aforementioned equation 
becomes: 



.56) 



Thus, the acceleration condition for a universe with quintom DE in the bulk reads: 



PB<-(A- PB ) or <j) 2 -& 2 < A/2 + V((f),cr). 



.57) 



Now we consider the 55-component of Einstein equations at the position of the brane, which leads to 
the Raychaudhuri equation 



*3 

6 



36 L 



5 b [3p b - pb) + p b ( 3p b + Pb 



(8.58) 



Using (|8.50p one can rewrite the above equation as 



a k 
a 6 



r '<A + ^5 2 
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k 4 

7^$b(Pb + 3p b 



i.59) 



Comparing equation (|8.55p with (|8.59p we provide a constraint on the bulk energy momentum tensor: 



(Sp B -Pb)= 4T|, 



(8.60) 



which for the quintom field with the energy-momentum tensor (|8.48p leads to 



(j) 2 — (J 2 = 0. (8.61) 

Therefore, we deduce that the time-dependent bulk quintom field influences the brane as a time- 
dependent cosmological constant, which can be written in terms of the quintom potential: ps = ~Pb = 
—V{4>) [H]. Thus, in order to acquire an accelerating universe the potential V must be a positive func- 
tion of time, and in order to be compatible with observations it should be decreasing with time. For 
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the particular solution of (|8.6ip in which eft and a are constant on the brane, we arrive at the natural 
cosmological constant on the brane, which in this case is induced by time-dependent bulk fields. 



Let us now consider a y-independent bulk quintom field with Lagrangian (j8.39[) [ 130|341j . The general- 
ized Friedmann equation in this case reads: 



1 + 



5 b K 



6p 2 



H 2 + - 
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where p? = SirG^ since v c — ^ \ • The brane Friedmann equation (j8.50p can be derived from fj8.62|) 
by letting go to infinity. We are interested in studying the effect of the quintom field on the brane. 
Ignoring the matter on the brane and the cosmological constant A, (|8.62f) can straightforwardly be 
rewritten as: 




.63) 



where ps is the energy density of the bulk quintom field on the brane, given by the first relation of 
(j8.48j) . The two different possible e-values (e = ±1), correspond to two different embeddings of the brane 
into the bulk spacetime [342 343J. Since the bulk quintom field satisfies the usual energy momentum 
conservation law on the brane (|8.54p . we have ps °c a~ 3 *- w+1 ^. Integrating (|8.63p for k = and w > — 1 
(that is (j) > &), shows that the scale factor a diverges at late times. Thus, the energy density of the 
bulk matter goes to zero for late times and thus it reaches a regime where it is small in comparison 
with 1/r 2 [341] . In the case w < — 1 (that is <j) < a), integrating (|8.63p indicates a vanishing scale factor 
a at late times, therefore the matter density goes to zero and we can use k\pb "C l/ r c- ^ n summary, 
in DGP model with a quintom DE fluid in the bulk space, one can expand (|8.63p under the condition 
that k^pb <C for the whole w-range. 

At zeroth order and for spatially flat metric, two different results, depending on the value of e, can be 
derived. Considering the case e = — 1 yields 



H 2 = 0, 

which describes an asymptotically static universe. Considering e = 1 leads to 



.64) 
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—z, or a(t) oc exp ( — 
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This branch provides the self acceleration solution in the late universe and it is the most important 
aspect of the model at hand. Therefore, the late-time behavior of the universe does not change even if 
we ignore the matter field on the brane and consider a scenario with the bulk quintom fields only. 

Let us now consider the "opposite" case, and neglect the bulk matter, considering only quintom fields 
confined on the brane, in DGP framework, with a Lagrangian as 



r<mat 



(8.66) 
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The energy-momentum tensor of the quintom fields on the brane is given by 



T, 



(8.67) 



In absence of bulk field and brane tension the first integral of the 00 component of (|8.67|) leads to 



in which is quintom energy density obtained by (|8.67p . Since the energy-momentum tensor on the 
brane is conserved, we could use the aforementioned procedure to examine the late-time cosmology on 
the brane. Integrating (|8.68p for a flat geometry indicates that for w > — 1 the scale factor a diverges at 
late times [309], while for w < — 1 the scale factor a vanishes at late times [130.341J. Thus, the energy 
density of quintom DE goes to zero for late times and reaches a regime where it is small in comparison 
with \jr 2 c . Expanding equation (|8.68|) under the condition fj?pb <C 1/r 2 provides an asymptotically 
static universe, H = 0, in the case e = —1 and a self accelerated phase, H = — , in the case e = 1. 

In summary, the presence of quintom DE on the brane or in the bulk does not change the late time 
behavior of the universe. In both cases and for the whole range of w (w < — 1 or w > —1), one can 
derive the self accelerating universe at late times. 



8.3 Other modified- gravity models 



We close this section mentioning that a quintom scenario can also be obtained in many other modified 
gravity models, namely, gravitational models with Gauss-Bonnet corrections [89 345 346 347 348 34 9~|350] . 
f{R) models with singular w = — 1 crossing [351) 352 353J, scalar-tensor models with nonminimal grav- 
itational couplings [128|354)355|356j . gravitational ghost condenstate models [200 357 358 359 360J (see 
Ref. [361J for a review). The analysis of the perturbations of these models can be studied by directly 
expanding metric fluctuations [362 363J or by the Parameterized Post-Friedmann approach [364 365. 366 J 
under certain cases, when one wants to check their rationalities. In the framework of these modified 
gravity models, the current cosmic acceleration may be achieved naturally. 



9 Energy conditions and quintom cosmology in the early universe 



In this section we investigate the violations of energy conditions in quintom cosmology, and we study 
its implications in the early universe. 



[309l344j 




(8.68) 
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9.1 Null Energy Condition 



It is well known that energy conditions play an important role in classical theory of general relativity and 
thermodynamics [367]. In classical general relativity it is usually convenient and efficient to restrict a 
physical system to satisfy one or more of energy conditions, for example in the proof of Hawking-Penrose 
singularity theorem [368 369 , the positive mass theorem |370|371] etc. Furthermore, in thermodynam- 
ics the energy conditions are the bases for obtaining entropy bounds [372)373) . Among those energy 
conditions, the null energy condition (NEC) is the weakest one, and it states that for any null vector 
n^ the stress-energy tensor T au should satisfy the relation 



Usually, the violation of NEC may lead to the breakdown of causality in general relativity and the 
violation of the second law of thermodynamics [374J. These pathologies suggest that the total stress 
tensor in a physical spacetime manifold needs to obey the NEC. In the framework of standard 4D FRW 
cosmology the NEC implies p + p > 0, which in turn gives rise to the constraint w u > — 1 on the EoS 
of the universe w u , defined as the ratio of pressure to energy density. 

In the epochs of universe evolution in which radiation is dominant the EoS of the universe w u is 
approximately equal to 1/3, while in the matter-dominated period w u is nearly zero, thus NEC is 
always satisfied. However, when the DE component is not negligible, with the NEC being satisfied, we 
acquire 

W u = Wm^AI + Wde^DE > ~1 , (9-2) 

where the subscripts 'M' and l DE' stand for matter and DE, respectively. With wm = 0, inequality 
(j9.2[) becomes 



From this inequality we can deduce that it is impossible for a universe to satisfy NEC if one of its 
component does not. In this case, a quintom scenario could be obtained if NEC is violated for a short 
period of time during the evolution of the universe. 

9.2 Quintom Bounce 

A bouncing universe with an initial contraction to a non-vanishing minimal radius and a subsequent 
expanding phase, provides a possible solution to the singularity problem of the standard Big Bang 
cosmology. For a successful bounce, it can be shown that within the framework of standard 4D FRW 
cosmology with Einstein gravity, the NEC is violated for a period of time around the bouncing point. 
Moreover, for the universe entering into the hot Big Bang era after the bouncing, the EoS must transit 
from w < — 1 to w > —1. 

We start with an examination on the necessary conditions required for a successful bounce |375|376|377|378)379)380|3"8T1 



T uv n"n u > . 



(9.1) 



WDE^DE > -1 • 



(9.3) 
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(we refer to [395J for a review on bounce cosmology). During the contracting phase the scale factor a(t) 
decreases (d < 0), while in the expanding phase it increases (d > 0). At the bouncing point, d = 0, 
and around this point a > for a period of time. Equivalently, in the bouncing cosmology the hubble 
parameter H changes from H < to H > 0, with H = at the bouncing point. A successful bounce 
requires 



*=-2M|(' ,+ P )=-2Mf( 1 + ! ") >0 - < 9 - 4 ' 

around this point and thus uu < — 1 in its neighborhood. After the bounce the universe needs to enter 
into the hot Big Bang era, otherwise it will reach the Big Rip singularity similarly to the phantom DE 
scenario [18] . This requirement leads the EoS to transit from w < —1 to w > —1, which is exactly a 
quintom scenario [391J. 



9.2.1 A phenomenological analysis 

Let us examine the possibility of obtaining the bouncing solution in a phenomenological quintom models, 
in which the EoS is described by: 



In (|9.5p r and s are parameters and we require that r < 1 and s > 0. One can see from (|9.5j) that w runs 

from negative infinity at t = to the cosmological constant boundary at t = y^x^f an d then it crosses 

this boundary. Assuming that the universe is dominated by the matter with the EoS given by (|9.5p . 
we solve the Friedmann equation, obtaining the corresponding evolution of hubble parameter H (t) and 
scale factor a(t) as: 



t 



a(t) 



(1 - r)t 2 + s 
i 



t 2 + 



1 



3(1-1") 



(9.6) 
(9.7) 



Here we choose t = as the bouncing point and we normalize a = 1 at this point. Thus, our solution 
provides a universe evolution with contracting (for t < 0), bouncing (at t = 0) and expanding (for 
t > 0) phases. In Fig. [30] we depict the evolution of the EoS, the Hubble parameter and the scale 
factor, as it arises from numerical elaboration. As we observe, a non-singular bounce is realized at t = 0, 
at a minimal non-vanishing scale factor o, with the Hubble parameter H running across zero. At the 
bouncing point w approaches negative infinity. 



9.2.2 Double-field quintom model 

Having presented the bouncing solution with the phenomenological quintom matter, we now study the 
bounce in a two-field quintom model with the action given by 
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Fig. 30. (Color online) Evolution of the EoS w, the Hubble parameter H and the scale factor a as a 
function of the cosmic time t. For the numerical elaboration we take r = 0.6 and s = 1. From Ref. $391f . 

^d^d^cj) - -d^o-d^a - V{</>, a) , 

where <j> is the canonical and a the phantom field. In the framework of FRW cosmology as usual we 
have p = ^(j) 2 — r^a 2 + V and p =\§ 2 — \& 2 — V and the cosmological equations read 

-\o- 2 + v) (9.8) 
(9.9) 
. (9.10) 

Thus, we see that a 2 = <fi 2 + 2V when H crosses zero, and moreover from (j9.4|) we deduce that a 2 = 4> 2 
when w crosses —1. These constraints can be easily satisfied in the parameter space of the model at 
hand. 

An interesting scenario realized by this model is a universe originated from a contraction, that has 
undergone a smooth bounce, then an inflationary period, and finally it returns to the standard thermal 
history. Models of this class have been studied in Refs. [385 391 396 397] and the dynamics of their per- 
turbations have been analyzed in Refs. [398 399 400 4 011402) . In the following we will provide explicitly 
one example to illustrate this scenario. 

We consider a double-field model in which the potential is only a function of the field eft and of Coleman- 
Weinberg form [403J: 
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which takes its maximum value Af 4 /16 at = and vanishes at the minima where 4> = ±v. Therefore, 
the scalar field a affects the evolution only around the bounce but it decreases quickly away from it. 



In order to discuss the perturbations explicitly, we first examine the evolution of the background. In this 
model a contracting universe can be driven to reach a minimal size during which the universe evolves 
like a matter-dominated one, and then a quasi-exponential expansion is following. The process to link 
the contraction and expansion is a smooth bounce, and the evolution of the hubble parameter can be 
treated as a linear function of the cosmic time approximately. 



For the background initial conditions we assume that (f) lies at one vacuum (for instance —v) when the 
universe is contracting, and moreover that a is small enough in order to be neglected. In this case the 
field (j> oscillates around — v leading the universe EoS to oscillate around w = 0, and so this state is in 
average similar to a matter-dominated one. Thus, the useful expressions of the background evolution 
write 



a~(- V ) 2 , H = -, H-rT 3 , (9.12) 
V 

where % = a' /a is the comoving hubble parameter and the prime denotes the derivative with respect to 
the comoving time 77. Moreover, since in the contracting phase the universe is dominated by the regular 
field and then we have an approximation 4>' 2 ~ 2(% 2 — %'), we can obtain another useful relation 



f. - 2m ~ U " (913) 

<// 2{n 2 - w) ' [ ' 

which will be used to calculate the metric perturbations. 



When the universe is contracting, the amplitude of 0-oscillations increases, while the contribution of 
the a- field grows rapidly. When the </>-field reaches the plateau, the bounce starts at the moment 
During the bounce we use the parametrization H(t) = a(t — ig) around the bounce point ts, and the 
coefficient a is a positive constant determined by numerical elaboration. This parametrization leads to 
a scale factor of the form a = \Z^~^rn - I n the bouncing phase the kinetic term of a reaches the 

maximal value, and from the equation of motion we deduce that (f)/(f) = —3H& 2 /(a 2 — -^q) — —3H 
when a is not very large. Finally, we obtain the approximate relations 

H - \kn - Vb), <f>" ^ -ZH<f/, |0| ~ e -if(™) 2 , (9.14) 
where we have defined y = 8aa 2 B /Ti. 



After the bounce, as the field 4> moves forward slowly along the plateau, the universe enters into an 
expanding phase at the moment ts+, and its EoS is approximately —1. Thus, the universe expands 
with its scale factor growing almost exponentially. In this phase we have the well-known relations for 
the background evolution 
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Fig. 31. (Color online) The evolution of (f)-field, of the EoS w and of the hubble pa- 
rameter H, in the model described in the text. For the numerical elaboration we choose 
the parameters A = 8.0 x 10~ 14 , v = 0.16M p; and we impose the initial conditions 
4> = -0.16M P , <p = 1.2 x Kr n M2, a = -0.14M P , & = 1.99 x l(T 14 M p 2 . From Ref. XSWj . 

a ~ , H ~ constant. (9.15) 

V 

Finally, when the </>-field reaches the vacuum state +v, it will oscillate again and the EoS of the universe 
will oscillate around zero as it happens before the bounce. 

In order to present this scenario explicitly, we perform a numerical elaboration and the results are 
depicted in Fig. [3TJ 

Finally, we mention that if the potential of the 4>-&eld is not flat enough, the inflationary stage could 
be very short. An extreme case of this scenario is that, after a matter-like contraction, the universe 
undergoes a smooth bounce and enters a normal expanding phase directly. This is the so-called "Matter 
Bounce" scenario, which has been intensively studied in the literature [404 405 40 6|407|408| . 

9.3 A cyclic scenario and oscillating universe 

The idea of cyclic universe was initially introduced in 1930's by Richard Tolman |409j . Since then there 
have been various proposals in the literature. The authors of Refs. |410|411|412|413j introduced a cyclic 
model in high dimensional string theory with an infinite and flat universe. Within a modified Friedmann 
equation the cyclic evolution of the universe can also be realized |414|415|416|417|418] . In Ref. [419] . it 
is shown that in the framework of loop quantum cosmology (LQC) a cyclic universe can be obtained 




77 



with the quintom matter. In this section, however, we will study the solution of oscillating universe in 
the absence of modifications of the standard 4D Einstein Gravity, within a flat universe. 



To begin with, let us examine in detail the conditions required for an oscillating solution. The basic 
picture for the evolution of the cyclic universe can be shown below: 

. expanding . contracting . 

...bounce > turn-around > bounce... . 

In the 4D FRW framework the Einstein equations can be written as: 

H^(^Y = -^ and ^ = -^#, (9-16) 
\aj 3M 2 a 6M 2 V ' 

where we have defined M 2 = g^j- H stands for the Hubble parameter, while p and p represent the 
energy density and pressure of the universe respectively. By definition, for a pivot (bounce or turn- 
around) process to occur, one must require that at the pivot point d = and a > around the 
bouncing point, while a < around the turn-around point. According to ()9.16f) . one obtains 

p = , p < (or p > 0) for the bounce (or turn-around), (9-17) 

or equivalently w = — > — oo (or + oo) at the bounce (or turn-around) point, with the parameter 
w being the EoS of the content of the universe. Thus, when the universe undergoes from bounce to 
turn-around w evolves from — oo to +oo, while in the reverse case w goes from +oo to — oo. This 
behavior requires w to cross over the cosmological constant boundary (w = —1) in these processes, 
which interestingly implies the necessity of the quintom matter in order to achieve the realization of the 
oscillating universe in 4D Einstein Gravity. 

An additional interesting feature of an oscillating universe is that it undergoes accelerations periodically, 
avoiding the Big Rip and Big Crunch, and furthermore we are able to unify inflation and current 
acceleration. We mention that the scale factor keeps increasing from one period to another and thus 
we are naturally led to a highly flat universe. This scenario was first proposed in Ref. [80], in which a 
parameterized quintom model was used and the coincidence problem was argued to be reconciled. 



9.3.1 A solution of oscillating universe in the double- field quintom model 



We consider the simplest quintom model consisting of two scalars with one being quintessence-like and 
another the phantom-like, with the usual action 



S 



d XyJ—g 



(9.18) 



As mentioned in previous sections the energy density and pressure are p = ^(ft — + V(4>,a) and 



p = — !<r 2 — V((j), a), and the equations of motion for these two fields write: 4> + 3H<f> + = and 
a + 3H& -V a = 0. 



Phenomenologically, a general potential form for a renormalizable model includes operators with di- 
mension 4 or less, including various powers of the scalar fields. We impose a Z2 symmetry, that is the 
potential remains invariant under the simultaneous transformations cp —> —(ft and a —> —a. Thus, the 
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potential of the model writes: 

V((ft, a) = V + \^n\(ft 2 + ^m 2 a 2 + 71 ^ + 72a 4 + gi (fta + g 2 (fta 3 + 53</>V + g^a. (9.19) 

From condition (|9.17p we can see that at both the bounce and the turn-around point a 2 = (ft 2 + 2V. 
Similarly, we deduce that p = —2V. When the universe undergoes a bounce the pressure is required to 
be negative, which implies the potential to be positive. However, when a turn-around takes place, the 
pressure of the universe is required to be positive, and consequently the potential must be negative. 
Therefore, the potential must contain a negative term in order to give rise to an oscillating scenario. 

In the model at hand the two scalar fields dominate the universe alternately, since the evolution from 
the bounce to the turn-around requires the transition from the phantom-dominated phase into the 
quintessence-dominated one, or vice versa. However, as pointed out in |81|82| this process does not 
happen if the two fields are decoupled, and thus interaction between the two fields is crucial. For a 
detailed quantitative study we consider 

V {(ft, a) = (A + \<fta) 2 + X -m 2 <ft 2 - ~mV . 

This potential acquires a negative value when (ft is near the origin and a large. However, due to the 
interaction the potential is still bounded from below and it is positive definite when the fields are both 
away from zero. This potential allows for the analytic solution 

(ft = \J~Aq cos mt , a = \J~Aq sin rat (9.20) 

with A = 2jPp , and with the parameter Aq describing the oscillation amplitude. In Fig. [32] we depict 
the evolution of the energy density and scale factor. 



9.3.2 Classifications of the solutions 

Let us study the detailed cosmological evolutions of the cyclic quintom model. First of all, from (|9.16|) 
one acquires the Hubble parameter 

H = ^ r (A + A 1 sin2mt) , (9.21) 

where Ai = ^^Aq. For different parameters, the evolution of the universe can be classified into five 
cases. Taking M p = 1 we have the following: 

Case (I): Ao = 0. In this case the Hubble parameter is given by H = sin2mt. Thus, the scale factor 
is 

lnaoccos2mt, (9.22) 

and therefore there is no spacetime singularity. In Fig. [33] we depict the evolution of the scale factor, 
the Hubble parameter, the energy density and the EoS. 

Case (II): < Ao < Ai. In this case the scale factor is 

In a oc C\t + C 2 cos 2mt , (9.23) 
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Fig. 32. (Color online) Evolution of the energy density and scale factor in a cyclic quintom universe. For 
each cycle the quintessence-like and phantom-like components dominate alternately. From Ref. \120j . 

where C\ = Aq/\^3M p , Ci = — Ao/8Mp. This solution also describes a cyclic universe, but both the 
minimal and maximal values of the scale factor increase cycle by cycle. Therefore, the average size of 
the universe is growing up gradually without Big Crunch or Big Rip singularities, although its scale 
factor experiences contractions and expansions alternately. Finally, the backward time evolution cannot 
lead to a shrinking scale factor either. These features are presented in Fig. [3H The graph for energy 
density and Hubble parameter is similar to FigJ33|. but the oscillating amplitude of the energy density 
in the expanding phase is larger than that os the contracting one. 

Case (III): Ao > Ai. The solution for the scale factor is the same as Case (II), however the evolution of 
the universe is different. The universe lies in the expanding period forever and there is no contracting 
phase, that is the accelerating expansion is periodical, without any bounce or turn-around. Since the 
EoS w is oscillating around "—1", the Big Rip can be avoided and the energy density and Hubble 
parameter are always positive. Furthermore, for reasonable parameters one can unify the inflationary 
period and the late-time acceleration. This behavior is shown in Fig. 

Case (IV): — |Ai| < Ao < 0. This case corresponds to a cyclic universe with decreasing minimal and 
maximal scale factor for each epoch. Contrary to Case (II), the total tendency is that the scale factor 
decreases cycle by cycle, since the contracting phase is longer than the expanding phase, but it never 
reaches zero. These features are depicted in Fig. [36] 

Case (V): Ao < — |Ai|. This case describes an enterally contracting universe, which corresponds to the 
reverse of Case (III), amd thus it is unphysical. 

Having presented the possible solutions, we make some final comments. One may notice that the EoS of 
the model oscillates around the cosmological constant boundary, and every time it approaches w = — 1 
it remains close to that for some time. Correspondingly, the Hubble parameter evolves near its maximal 
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Fig. 33. (Color online) Evolution for Case (I). This figure presents an exactly cyclic universe. The 
scale factor oscillates between the minimal and maximal value. In the numerical elaboration we take 
m = 3 x 10 -3 and Aq = 300 or equivalently A\ ~ 0.39. The Hubble parameter is depicted by the purple 
line, while the energy density by the red line. From Ref. [120]. 




Fig. 34. (Color online) Evolution for Case (II). m and Aq are the same as those used in Fig. [23 and 
A = 0.10. From Ref. \TWj . 

value and makes the scale factor expand exponentially. This period corresponds to an inflationary stage 
after the bounce. It is a significant process since it dilutes the relics created in the last cycle. Furthermore, 
entropy can be diluted too, providing a solution to the problem of enteral entropy increase in cyclic 
cosmology. Additionally, some of the primordial perturbations are able to exit the horizon and re-enter 
the horizon when the inflationary stage ceases. When these perturbations re-enter they lead to new 



81 





Fig. 36. (Color online) Evolution for Case (IV). The values of m and Aq are the same with Fia l33l and 
A = -0.10. From Ref. UWj . 

structures in the next cycle. 



Lastly, let us constrain the model parameters (m, Ao, ^4o)- First of all, we require the oscillation period of 
the Hubble parameter to be no less than 2 times the age of our universe which is of the order of the present 
Hubble time. Thus, from (19T2TJ) we obtain the period T = ^ ~ 0(Fo _1 ), where H^ 1 ~ lO^M" 1 , and 
therefore m ~ O(10~ 60 )M p . Secondly we require the maximal value of the Hubble parameter to be able 



S2 



to reach the inflationary energy scale. The maximal value of H is given by (Ao + Ai), when (sin2mi) 
reaches its maximum. For case (I), (III) and (IV), |A | < Ai, so we have (A + Ai) ~ 0( r jfy). If we 

consider inflation with energy density, for example around 0(1O~ 2O )M P 4 , we find that Aq must be of 
0(lO 5O )M p 2 . Such a large value of Aq indicates that the scalar fields 4> an d o take values above the 
Planck scale, which makes our effective lagrangian description invalid. One possibility of solving this 
problem is to consider a model with a large number of quintom fields. In this case the Hubble parameter 
is amplified by a pre-factor y/~N with N the number of quintom fields. If N is larger than 10 50 , Aq can 
be relaxed to 0(1)M 2 or less. 



10 Concluding remarks 

Since discovered in 1998, the nature of DE has become one of the most intriguing puzzles of modern 
physics and it has been widely investigated. The simplest candidate of DE is a cosmological constant, 
but it suffers from the well-known fine-tuning and coincidence problems. Alternatively, dynamical DE 
models have been proposed, such as quintessence, phantom, and k-essence. Since at present we know 
very little about the theoretical aspects of DE, the cosmological observations play a crucial role in 
enlightening our picture. In particular, since astronomical observations have shown a mild preference 
for an EoS of DE smaller than —1 at present, a phantom model appeared in the literature. 

In these lines, with accumulated observational data a scenario of quintom, with DE EoS crossing — 1 
during evolution, has been constructed in the literature. If such a class of dynamical DE scenario were 
verified by future observations, it would be a challenge to the model-building of DE. This is because of 
a No-Go theorem, which forbids dynamical models with a single scalar field to lead to an EoS crossing 
over the cosmological constant boundary in the frame of Einstein's gravity. The current work aims 
at presenting a review of successful examples of quintom models and the corresponding observational 
consequences. 

We studied in detail the simplest quintom model, which involves two scalar fields, with one being 
quintessence and the other being phantom. We have shown that these models possess stable critical 
points in the phase space, and in addition certain dualities appear, with form-invariant cosmological 
equations. 

As a consistent consideration we analyzed the behavior of quintom perturbations. In the example of the 
simplest quintom model there are two degrees of freedom in its perturbations, which are finite and stable. 
It should be noticed that since the quintom scenario usually involves extra degrees of freedom, additional 
relative pressure perturbations are inevitably produced in order to seed entropy perturbations. 

We also provided alternative approaches to realize the quintom scenario in a scalar field system within 
standard Einstein's gravity, in which higher derivative terms were introduced. Specifically, we con- 
structed a Lee- Wick and a string inspired quintom model. For the latter, the perturbation propagations 
are quite different than usual, with a time-dependent sound speed due to the non-canonical kinetic term. 

Usually, a quintom model constructed by a double-scalar system and a single scale with higher deriva- 
tives, suffers from the problem of ambiguous quantum behavior inherited from the DE phantom model. 
This feature offered a motivation to study the possibilities of realizing a quintom scenario without a 
ghost mode. An attempt on solving this problem is to construct a string-inspired quintom model, but 
it is still unclear how to be quantized since it involves highly non-linear terms. 
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Furthermore, we considered the model-building in the framework of a non-scalar system. A model of 
spinor quintom is able to evade the drawbacks of considering a phantom field, since its kinetic term 
is well defined, and the sign-change of its potential derivative with respect to the scalar bilinear ipip 
can realize the — 1-crossing. However, it is non-renormalizable if we consider higher order effects at the 
quantum level. 

Another interesting issue is the investigation of a quintom scenario in braneworld cosmology, which is 
motivated by higher dimensional constructions, closer to fundamental theories. Braneworld models were 
studied actively in recent years, offering a new path to describe DE. In this review we presented the 
possibility of realizing a quintom scenario in DGP braneworld model. In such a construction, solutions 
of late-time super-acceleration can be obtained when matter components are driven by a manifestation 
of the excruciatingly slow leakage of gravity into the extra dimension. In particular, we showed that 
for a canonical field living on the brane the —1-crossing takes place from above to below, while in the 
phantom case it takes place from below to above. 

Concerning the theoretical implications of quintom cosmology, an important issue is the violation of 
energy conditions. Such conditions usually play a crucial role in classical theory of general relativity. An 
important predication of general relativity is that a singularity cannot be avoided under assumptions 
of certain energy conditions, which was originally proved by Penrose and Hawking and later developed 
into cosmological framework by Borde, Vilenkin and Guth |420j . Therefore, in a quintom scenario with 
NEC violation, one may obtain a bouncing solution at early times of the universe. This may provide 
a possible solution to the singularity problem of standard Big Bang cosmology. Specifically, we studied 
bouncing cosmology realized by a model of double-field quintom. To extend, we also provided a scenario 
of oscillating universe, in which it experiences expansions and contractions periodically. 

As an end, we would like to comment on some unsettled issues in quintom cosmology. Although the 
quintom paradigm has been studied intensively in recent years, its nature is still unclear. In the following, 
we list three main questions, which are crucial to the developments of quintom cosmology. 

• (a) how to combine quintom cosmology with the well-known particle physics or fundamental theories? 

• (b) how to construct a ghost-free quintom scenario with robust quantum behavior? 

• (c) why the DE EoS crosses the cosmological constant boundary around today (the so-called second 
coincidence problem)? 

Recently, physicists have made many attempts on partially addressing the above questions in quintom 
cosmology. For example, a quintom scenario can be realized by a rolling tachyon in the cubic string field 
theory |101|421|422| due to a non-local effect |423j . Interestingly, motivated by studies on thermody- 
namics of a cosmological structure [424 425 426 427 428 429 430], a cosmic holographic bound (originally 
suggested by |431)432| and comprehensively reviewed in [433J) on dark energy dynamics allows a period 
of quintom scenario with NEC violation |434I435I436I437I438I439] . 

However, the three aforementioned questions remain open, requiring to be faced by the coming theoret- 
ical studies on quintom cosmology. In the meantime, cosmological observations of higher accuracy are 
needed in order to determine whether the DE EoS has crossed the cosmological constant boundary. If 
this is indeed true, cosmology will enter in a very interesting and challenging era. 
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